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Abstract: Black hole uniqueness is known to fail in higher dimensions, and the multi- 
plicity of black hole phases leads to phase transitions physics in General Relativity. The 
black-hole black-string transition is a prime realization of such a system and its phase dia- 
gram has been the subject of considerable study in the last few years. The most surprising 
results seem to be the appearance of critical dimensions where the qualitative behavior of 
the system changes, and a novel kind of topology change. Recently, a full phase diagram 
was determined numerically, confirming earlier predictions for a merger of the black-hole 
and black string phases and giving very strong evidence that the end-state of the Gregory- 
Laflamme instability is a black hole (in the dimension range 5 < D < 13). Here this 
progress is reviewed, illustrated with figures, put into a wider context, and the still open 
questions are listed. 
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To Dorit, Inbal and Neta, 
my wife and daughters 

1. Introduction 

In this introduction wc will present the general background for studying General Relativity 
in higher dimensions and the novel field of phase transitions in General Relativity. We will 
list the systems known to exhibit phase transitions, and take the opportunity to discuss the 
rotating black ring before we proceed to concentrate on our system of choice, the black-hole 
black-string transition. 

Why study GR in higher dimensions? There are several good reasons to study 
General Relativity (GR) in higher dimensions, namely D > 4, where D is the total space- 
time dimension. From a theoretical point of view there is nothing in GR that restricts us 
to D = 4. On the contrary, the theory is independent of D, and D should be considered 
as a parameter. It is common practice in theoretical physics to explore large regions of 
parameter space of a theory in order to enhance its understanding, rather than restrict 
to the experimental values and GR should be no exception. For example, in the study 
of gauge theories it is standard to consider various possibilities for the gauge group and 
matter content which differ from the standard model. 

Additional reasons to study higher dimensional GR include string theory and the phe- 
nomenological scenario of "large extra dimensions", as we proceed to discuss. String theory 
has a "built-in" preference for higher dimensional spacetimes with 10 (the "critical dimen- 
sion") or 11 dimensions, where the extra dimensions must be compactified. This preference 
originates in the cancellation of the conformal (quantum) anomaly in lOd which is neces- 
sary for the consistency of weakly coupled string theories. The "large extra dimensions" 
scenario (which is presumably string theory inspired) stresses the following important re- 
alizations: that to date gravity is measured only down to l/x-lmm range (which is an 
"astronomically" poor resolution relative to the one we have for other forces), that it is 
quite consistent to assume the existence of a compact dimension(s) smaller than the ex- 
perimental bound and that the situation can be rectified only by improving gravitational 
and accelerator experiments. 

The novel feature - non-uniqueness of black objects. Often when we generalize a 

problem to allow for an arbitrary dimension the qualitative features do not change and thus 
the generalization does not produce "new physics" , even if the quantitative expressions are 
different. However, in GR we do find qualitative changes. If we roughly divide the field of 
General Relativity into black holes, gravitational waves and cosmology, we find a qualitative 
change in the first of these categories: one of the basic properties of 4d black holes changes, 
namely black hole uniqueness.^ 

^ Other qualitative differences include the disappearance of stable circular orbits for D > 4 (in Newtonian 

gravity), the absence of propagating gravitational waves in D < 4, and the Belinskii-Khalatnikov-Lifshitz 
(BKL) analysis of the approach to a space-like singularity, where there is a critical dimension Dgj^j^ = 10, 
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Here we should digress to make the distinction between two closely related black hole 
notions: "no hair" and "uniqueness" (see for example "No hair" denotes the feature 
that the space of black hole solutions has a small dimension usually parameterized by 
asympototically measured quantities (mass, angular momentum and electric charge, for 
example) much like macroscopic thermodynamical variables. In this respect a black hole 
strongly contrasts with a non-black-hole star which typically has a much larger number of 
characteristics such as its internal matter ingredients each with its own equation of state and 
spatial distribution possibly resulting in an unbounded number of independent multipoles 
for mass, charge and angular momentum. Whether this "no-hair" property continues to 
hold for higher dimensional black holes could depend on the way one chooses to generalize 
it. If one generalizes "no hair" to mean that the solutions are determined in term of a 
small number of (not necessarily conserved) asymptotic data then it continues to hold in 
higher dimensions as far as we know. However, if one would choose the more restrictive 
definition which requires conserved charges then this property fails in higher dimensions as 
was demonstrated in the generalized rotating black ring ||3|, ^, ^ , whose parameters include 
some non-conserved dipole charges. 

"Uniqueness" on the other hand is the more specialized statement that a choice of all 
of these asymptotic black hole parameters selects a unique black hole rather than a discrete 
set. In other words, that only a single branch of solutions exists. In 4d uniqueness was 
proven to hold, namely that given the mass, charge and angular momentum (satisfying 
some inequalities to ensure the existence of a solution with no naked singularities) there 
is a unique black hole. However, the proof relies heavily on properties which are special 
to 4d: Hawking's proof that the horizon topology has to be and the simplifying gauge 
choices of Weyl-Papapetrou and Ernst (see Q for references to original papers and reviews 
and for a speculative generalization of uniqueness to higher dimensions). See ||^, ^ for a 
determination of the allowed horizon topologies in certain D > A. 

The breakdown of black hole uniqueness in higher dimensions implies the coexistence 
of several phases with the same asymptotic charges on a non-trivial phase diagram. Phase 
transitions between the various phases should occur as parameters are changed. As always 
one may define the order of the phase transition. It could be a first order transition in which 
case it is triggered non-perturbatively by a competition of entropies between two phases 
which are separated by a finite distance in configuration space, or it could be of second or 
higher order, in which case it is triggered by perturbative tachyons and the transition is 
smooth (see subsection 3.2). 

Such first order transitions would be accompanied by an exceptional release of energy, 
sometimes called a thunderbolt,'^ simply since the total mass of the final state must be 
lower than or equal to that of the initial state and the excess energy must be lost through 
radiation. Moreover, exact mass equality is highly unlikely, but rather a loss of mass is 
natural as spacetime would undergo violent changes including sometimes the roll-down of 
a tachyonic mode. 



such that for D > Dgj^j^ the system becomes non-chaotic (see the review |l| and references therein). 

^This term was introduced by |^] for a certain gravitational shock wave in the presence of a naked 
singularity and seems appropriate for the system under study as well. 
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Figure 1: In asymptotically flat 5d spacetime uniqueness is violated by the co-existence of the 
rotating black hole and the rotating black ring. The figure shows the range of existence of each 
phase on the dimensionless angular momentum axis. Note that for 1 < nJ'^/M^ < 32/27 three 
phases co-exist. 

The two systems. To date we know of two systems with higher dimensional non- 
uniqueness resulting in non-trivial phase transition physics 

• The rotating ring 

• The black-hole black-string transition 

The latter was chosen for a thorough study of its phase structure which is the subject 
of this review, presumably since it is somewhat simpler to analyze on account of the smaller 
number of metric functions and its higher degree of symmetry. Before proceeding to analyze 
it in detail, we discuss the other example, the black ring. 

The rotating black ring lives in the flat (and topologically trivial) 5d background. 
Spherical rotating black holes solutions in higher dimensions, which generalize the 4d Kerr 
solution were already found in 1986 by Myers and Perry These solutions have an 
horizon topology in 5d and display a maximal angular momentum (at fixed mass). In 2001 
Emparan and Reall discovered a beautiful solution, the ring, with horizon topology 
X S^, and with an angular momentum which is bounded from below, but not from 
above (at fixed mass). Figure |l] shows the regions on the angular momentum axis which 
are occupied by the various phases, and one notes that there is a middle region where 
three phases coexist - one black hole and two rings. Unlike the black string, rotating 
ring solutions are known only in 5d, presumably due to the special property that in 5d 
the centrifugal potential and the (Newtonian) gravitation potential have the same (1/r^) 
r-dependence. 

For some time it was not clear whether the black ring is stable, and despite some recent 
findings the issue is not settled yet. In 2004 charged rings were shown to be BPS [^2| and 
hence plausibly "super-stable" (that is, non-perturbatively stable). Although the stability 
of the original non-BPS ring of is still undetermined (however, see |13] for interesting 



partial results on stability), there is comfort in knowing that some of its closest relatives 
which share many of their outstanding properties are plausibly stable. Soon after, several 
groups made progress in obtaining larger families of ring solutions, both BPS H, |, |, |l5l 
and non-BPS ||l^, all of them restricted to 5d. Moreover, the inclusion of non-conserved 
dipole moments in [|| demonstrates the "no-hair" principle in higher dimensions must be 
generalized at least to allow for non-conserved quantities. 
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As the recent discovery of families of black rings demonstrates, the black ring may 
hold further surprises. In particular, we do not know the full parameter space for rings, 
and we know close to nothing about the associated phase transitions. Thus rotating rings 
constitute a promising and active field of research. 



Outline. At this point we set aside the topic of black rings until the discussion section and 
we turn in the next section to the other example for non-uniqueness, the black-hole black- 
string transition which is the main topic of this review. In section ^ the physical set-up is 
described and the questions of interest are formulated. In section ^ we describe the analytic 
considerations that culminate in subsection 3^ to a certain suggestive qualitative form of 
the phase diagram which is compared there with numerical data. Section ^ describes the 
quantitative tools that were employed in order to obtain solutions, including both numerical 
and analytic methods. Finally, related work is described in section ^ and we conclude with 
a summary of the results and a discussion of open questions in section 0. 



2. Set-up and formulation of questions 
2.1 Background metric and phases 

Background metric. We consider a background with extra compact dimensions. In such 
a background one expects to find several phases of black objects depending on the relative 
size of the object and the relevant length scales in the compact dimensions. For simplicity 
we discuss here pure GR (the only field is the metric) with no cosmological constant. Thus 
the backgrounds considered are of the form M*^"^'^ x where is any p-dimensional 
compact Ricci-flat manifold, d is the number of extended spacetime dimensions, and the 
total spacetime dimension is D = d + p. 

The simplest compactifying manifold is a single compact dimension X = S^, and 
accordingly that was the X considered in most of the research so far. was chosen 
not only for its simplicity but also since while more involved X will have several phases of 
black objects, the phase transition physics between any two specific phases is expected to be 
essentially similar (generically) to the case. Some research was devoted to X = T^, the 
p-dimensional torus and we shall discuss it later. Other possibilities for X include 

K3 and Calabi-Yau threefolds, as well as Ricci-flat spaces X which are not supersymmetric. 

Thus we consider a background with a single compact dimension of size L, namely 
>^ gi ^ ^-j^fi D — ^_|_]^ > 5 ^-^T^Q lower bound on D is set in order to avoid spacetimes with 
2 or less extended spatial dimensions where the presence of a massive source is inconsistent 



with asymptotic flatness, see [19, EO, 21| for a limited analogue in D = 4). 



Black objects. The non-rotating black objects in which we are interested are static and 
spherically symmetric. Thus, the essential geometry is 2d after suppressing time and the 
angular coordinates in the extended dimensions. Our coordinates are defined in figure ^. 

Such solutions are characterized by 3 dimensionful parameters M, L and Gn where 
M is the mass of the black object (measured in the asymptotic d dimensional spacetime. 
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Figure 2: Definition of coordinates. For backgrounds with a single compact dimension tlie essential 
geometry is 2d after suppressing time and angular coordinates in the extended dimensions. The 
cylindrical coordinates (r, z) are defined such that z ^ z + L is the coordinate along the compact 
dimension and r is the radial coordinate in the extended spatial directions. For black holes we 
define another set of local coordinates {p, x), defined only for p < L/2, which are radial coordinates 
in the 2d plane with origin at the center of the BH. 



and the detailed expression is given in and Gn is Newton's constant. These define a 
single dimensionless parameter ^ 

GnM 

^i:=jj^. (2.1) 

Alternatively one may use a different parameterization of the solutions such as replacing 
M by P, the inverse temperature.^ Correspondingly one may define another dimensionless 
parameter 

/i^ oc ^ , (2.2) 
where the proportionality factor may be chosen later by convenience. In thermodynamic 



terms, the parameter (2.1) or (2.2) is the "control parameter" of the system, and the 
choice between these two depends on whether one prefers the micro-canonical or canonical 
ensembles, respectively. 

In this background one expects at least two phases of black object solutions: when 
/X <C 1, namely the size of the black object is small (compared to L the size of the extra 
dimension) one expects the region near the object to closely resemble a D-dimensional 
black hole, while as one increases the mass one expects that at some point the black hole 
will no longer fit in the compact dimension and a black string, whose horizon winds around 
the compact dimension will be formed. The precise distinction between these two phases 
is give by 

Definition : We distinguish between the black hole (BH) and the black string according 
to their horizon topology which is either spherical — S^~^ or cylindrical — S"^"^ x X, 
respectively. 

These phases are illustrated in figures |3|,^,|5|(b). We shall sometimes refer to such a black 
hole localized in a compact dimension as a "caged black hole". 



^We are dealing with classical GR, and thus we do not set h — 1. 

*More precisely in order to avoid using h we define here (3 = 2tv/k, namely the period of the Euclidean 
time direction, and k is the surface gravity. 
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Figure 3: The uniform black string, rg is its Schwarzschild radius. 



Applications. Before proceeding to discuss the phases in more detail, let us mention some 
applications that contribute to its importance, beyond its considerable intrinsic value. In 
String Theory it has attracted continued interest, particularly regarding the thermody- 
namic phase diagram for various gravity theories and/or field theories |2^, |2^] which are 
related by dualities to the higher dimensional origin of brane solitons (such as the M-theory 
origin of string branes), where the physics is significantly affected by the question whether 
they are localized in the compact dimensions or wrap them. Another field of application 
is black holes on brane- worlds |24, a problem closely related to the one discussed here, 
only the background in which the black objects live includes not only an extra dimension 
but also a "phenomenological" brane localized in that dimension and carrying the fields of 
the standard model. 



We now proceed to discuss the two phases with more detail. 
The black string. We can readily write down solutions which describe uniform black 
strings (see figure 

= ^^Schw + (2.3) 
where ds^ is the metric on X, which for our central example, an parameterized by the 
coordinate z, is just 



ds 



X 



dz^ , 



(2.4) 



and dsa i is the d-dimensional Schwarzschild black hole (also known as Schwarzschild- 
Schw ^ 

Tangherlini ||2^), which is given by 



'^^Schw 



-f{p)dt' + j^^dp' + p^dnj_. 



where 



fip) = 1 



41 



(2.5) 
(2.6) 



dJl? 2 is the metric on the sphere S' 



d-2 



dQl_2 = dx^ + sm^xd9j + + (sin^ x sin^ 6*1... sin^ 6d-4)d9j_2 



(2.7) 



Po is related to the black hole mass, M, via |T^] 

Po 



{d - 2) 



(2.8) 
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Figure 4: A caged black hole (BH). Newtonian equipotential lines are shown. 



where Gd is the d-dimensional Newton constant, and Vt^- 



a unit sphere S'^~^. The relation between ro and the inverse temperature /? is 



27r 



47r 



(2.9) 



/'(ro) d-3 

These metrics are Ricci flat as a result of being a direct product of Ricci flat metrics. 
They are called "uniform" for being a direct product with X (moreover, for the full 
metric is z-independent). Later we will encounter also non-uniform strings (see flgure 
|5|(b)). Note that for general p > 1 (namely dim{X)) these metrics actually describe p- 
branes rather than strings. 

The uniform black string solution is valid for any tq (and fixed X). However, we shall 
soon see that for "thin" enough strings, namely small enough tq, an instability develops. 



Caged black holes. One expects localized black hole (BH) solutions to exist (see figure 
Q), intuitively obtained by constructing a black hole locally without ever being "aware" of 
the compactness of the some of the dimensions, at least as long as the black hole is much 
smaller than the compact dimensions (and the number of extended spacetime dimension 
is d > 4 to avoid problems with asymptotics) . 

As the black hole grows it will start feeling the presence of the compact dimensions 
and it will deform accordingly. At some critical /i one may expect that the black hole will 
be too large to fit into X, and so the mass of this phase will be bounded from above. 

Unlike the uniform black string there is no explicit metric that we can write down. 



This situation was confronted by two methods: an analytic perturbative expansion [27, 28 



|29| , |30[ | and numerical analysis 33, 34 1. Both techniques will be described in section 

^ and here we only note that the analytic method is useful for small black holes (actually 
fif^ is the small parameter for the perturbation series), while for large black holes, where 
the interesting phase transition physics occurs the numerical methods are essential. The 
existence of both techniques created a healthy feedback where both methods were used to 
test and improve each other. 
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(a) (b) 

Figure 5: (a) The Gregory-Laflamme instability, (b) A non-uniform black string. 



2.2 Gregory-Laflamme instability 



Gregory and Laflamme (GL, 1993 |35]) discovered that the uniform black string solution 
( |2.5| ) develops a z-dependent metric-instability below a certain critical mass (see fig- 
ure ^). By a "metric-instability" one means that when one analyzes the spectrum of 
frequencies-squared for small perturbations around this background, a negative eigenvalue 
is found. ^ 

In hindsight, this instability makes a lot of sense. In general, gravity has a tendency to 
clump matter. For example, a uniform distribution of gravitating matter ("gas") is know 
to be unstable against the formation of inhomogeneities (the so-called "Jeans instability"): 
when an inhomogeneity forms the denser regions exert a stronger gravitational pull on 
their neighborhood, thereby triggering an unstable positive feedback. Similarly here, a 
long enough string "wants" to develop inhomogeneities (if it is short enough then it gets 
stabilized by the energetic costs of spatial gradients). Another perspective is to recall 
that the Schwarzschild black hole has negative specific heat (black hole thermodynamics). 
While this is not enough to de-stabilize a single black hole, it should certainly destabilize 
a homogeneous collection of black holes, namely a black string, which could increase its 
entropy by re-distributing its mass non-uniformly. This intuition is the basis for the Gubser- 
Mitra Correlated Stability Conjecture ||36| , |37| which states that a homogeneous black brane 
is (classically) perturbatively unstable if and only if the dimensionally reduced black hole is 
thermodynamically unstable (semi-classically)^. From now on we continue to discuss only 
perturbative instabilities. 

The main results of GL are summarized in figure |6| (taken from |]35| , |39[| ) which depicts 
the inverse decay times as a function of vq/L for total spacetime dimensions 5 < D < 10. 



^Such a metric instability is also known as a "tachyon", where the latter term is used in a more gen- 
eral sense than the "usual" 4d tachyonic field. While one usually considers 4d tachyonic fields cj>, whose 



Lagrangian behaves as ~ i \^'^ — (V^(^)^ -I- m^? (/!>^J , where m^? > and V stands for a 3d spatial gradient, 
one also generalizes it to arbitrary spatial dimension, including spatial dimension 0, which is the case here, 
when we take to be the amplitude of the GL mode and = Q? , where fi, the inverse decay time, is to 
be defined shortly in figure 6. 

®See \Rm for proofs of certain aspects of this conjecture. 
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Figure 6: Characteristic inverse decay times $7 (giving the perturbation an e^* time dependence) 
as a function of the perturbation wavenumber /i :— nro/L — rQk/2 (proportional to fip in our 
notation) for A < D < 9 where D is the extended space-time dimension (d in our notation) in 
backgrounds with an compactification. kcL, the critical Gregory-Laflamme wave-numbers, are 
the maximal wavenumbers for which the instability exists - namely, the intersection points with 
the horizontal axis. The bold points correspond to value calculated numerically and the lines have 
been traced to guide the eye (all lines converge to the origin = /i = which is a non-physical 
gauge mode). Reproduced from |9|. 



From figure ^ we see that the tachyonic mode appears for wavenumbers k (at fixed tq) 
which are lower than a critical wavenumber kcL (which depends on d). In order to find 
kchid) it is not necessary to look at perturbations in D dimensions, but rather it suffices 
to find the negative mode of the Euclidean d dimensional Schwarzschild black hole. This 
mode, discovered by Gross, Perry and Yaffe I^Oj (in the 4d case) and hence denoted here 
by hcpY satisfies 

Lschw hcPY = —^GPY hcPY (2-10) 
where Lgchw is the Lichnerowicz operator for perturbations in the Schwarzschild back- 
ground and —Xgpy is the negative eigenvalue. Given Hgpy tbe marginally tachyonic 
mode is given by |4l|, |38| ^ 

h = hcpY exp {i kcL z) 



kcL ■■= VXgpy (2.11) 

In [^] the critical GL lengths were obtained for Schwarzschild black holes in various 
dimensions (see table |^). Prom these the high d asymptotic form was extracted and later 
proven analytically in |l^ to be 

kcL^Vd—. (2.12) 

This means that for large d the black string becomes unstable at a compactification length 
Lgl = '^T^/kcL ~ To/^/d when it is quite "fat" (namely ro ^ Lgl) and indicates that such 
a string would not decay into a black hole which would not "fit" inside the extra dimension. 



^Footnote 4 of jssj explains "This relationship (eq. 2.11 - BK) was noted in ahhough the connection 
with classical stability was not appreciated." 
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At k = kcL the GL mode is marginally tachyonic, namely a zero-mode. Morse theory 
arguments strongly suggest^ that this zero-mode produces a branch of solutions emanating 
from the GL point describing non-uniform strings due to the z-dependence of the GL 
mode. 



d 


4 


5 


6 


7 


8 


9 


10 


11 


kcL 


.876 


1.27 


1.58 


1.85 


2.09 


2.30 


2.50 


2.69 


d 


12 


13 


14 


15 


19 


29 


49 


99 


kcL 


2.87 


3.03 


3.19 


3.34 


3.89 


5.06 


6.72 


9.75 



Table 1: Numerically computed static mode wavenumbers k^L in units of Tq as a function of d, 
the number of extended space-time dimensions [EtI 521 . 



The end-state. Whenever one discovers a perturbative tachyon indicating a decay, the 
question of its end-state is naturally raised. As the end-state configuration often lies away 
from the initial configuration a perturbative analysis does not suffice and one needs global 
information regarding all stable static solutions which is more difficult to obtain. 

Gregory and Laflamme believed the end-state to be the black hole, both since that 
was the only other phase they knew about and since comparing entropies at small /i one 
notices that the black hole phase has superior entropy in this regime. Indeed for small 
where the black hole is well approximated by a D dimensional spherical black hole, the 
entropies scale as 



Sbh ~ r-g^ ~ 



Sst ~ ~ , (2.13) 

in units where L = = G£) = 1. It is seen that the exponent {D — 2)/{D — 3) = 
1 + 1/(D — 3) is a monotonically decreasing function of its argument and hence the exponent 
is smaller for D > d (the black hole) resulting in a larger area. 



More recently Horowitz and Maeda |43] showed that the black string horizon cannot 



pinch in finite "horizon time" (namely, finite affine parameter along the horizon generators). 
They interpreted that as an indication that surprisingly a black hole could not be the end- 
state of decay and predicted instead the existence of a stable non-uniform string phase 
that would serve as an end-state. The argument for pinching in infinite horizon time 
relied on assuming the increasing area theorem for an event horizon and applying it to 
an area element at the "waist" - the inward collapsing region of the event horizon. The 
extension to the claim on the end-state involved estimates on why infinite horizon time 
should imply infinite asymptotic time (time for an asymptotic observer). While these 
claims stimulated much of the research reported here, and while numerical evidence lends 
support for "pinching in infinite horizon time" (see subsection ^]^) , strong evidence against 



*The phrase "strongly suggest" is used conservatively due to possible subtleties in the argument which 



are indicated in subsection 3.3 and were not explored in full rigor. However, it is the author's opinion that 



Morse theory arguments essentially guarantee the existence of the non-uniform branch. 
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the end-state being a non-uniform string will be described as this review proceeds, implying 
that the end-point is actually the black hole phase as originally argued by Gregory and 
Laflamme (at least in dimension D < 13). See the summary section for a more complete 
discussion. 

2.3 Issues 

Let us formulate some major issues or questions regarding this phase transition. These 
issues may be roughly divided into two groups: static and time evolution. 
The static issues include 

• End-state of decay. 

• Qualitative form of the phase diagram including the determination of all static phases. 

• Detailed quantitative data on the phase diagram: the domain of existence of each 
phase and the determination of critical points. 

During the last couple of years there was significant progress on the static issues, 
resulting also in the surprising discoveries of critical dimensions and a topology change. 
The deepest issues belong however to the time evolution 

• The spacetime structure, namely determination of the Penrose diagram, or an ap- 
propriate generalization thereof. 

• A naked singularity and a violation of Cosmic Censorship. 

It is plausible that as the black string pinches a naked singularity is formed, naively 
because the singularity which "originally" winds the compact dimension gets "bro- 
ken", perhaps at the event of pinching. Another argument comes from the clash 
between the arguments of Q| and results on the system's phase diagram 
Possibilities include a problem with the assumption that there are no singularities 
strictly outside the horizon and an infinite duration with respect to "horizon-time" 
(horizon affine parameter) while the asymptotic-time duration is finite. Note that 
the initial conditions in this case are generic, unlike known examples of naked singu- 
larities. 

• A thunderbolt and quantum gravity. 

The decay is accompanied by a release of energy (after all, a tachyon is involved) in 
the form of radiation (see ||4^). It is plausible that this radiation pulse is classically 
singular (a "thunderbolt"), perhaps due to its origin from the naked singularity. In 
such a case it is quite plausible that some knowledge of quantum gravity will be 
necessary in order to understand this outgoing radiation. 

While there was much progress on the static issues, there was practically none on the 
time evolution issues, and these remain unsolved. 



- 12 - 



3. Qualitative features 

In order to understand the phase transition physics and to resolve the issue of the end-state 
it suffices to map out ah static and stable solutions of the system, since the end-state is 
certainly static and stable. But actually, Morse theory arguments will lead us to consider all 
static solutions whether stable or not, in order to take advantage of a "phase conservation 
rule" which is a powerful qualitative constraint on the form of the phase diagram. So we 
seek the phase diagram of all static solutions as a function of /i, and throughout this review 
we will restrict ourselves to static aspects of the system.^ 

In this section we seek to determine the phase diagram qualitatively, and the quanti- 
tative aspects will be described in the next section. 



3.1 Order parameter 

We wish to define an order parameter A such that the uniform string will have A = and 
the emerging non-uniform branch (from the GL point) will have A 7^ 0, namely A should 
be a measure of non-uniformity^''. Actually, it is desirable to have both the black strings 
and the black holes at finite values of A, motivated by the expectation for a merger of the 



two due to Morse theory arguments as will be explained in subsection 3.3. 

It turns out that an asymptotic analysis of the metric and the associated charges 
furnishes physically meaningful candidates [46, However, it should be noted that the 



central discussion on the qualitative form of the phase diagram that will culminate in 
subsection ^ will be independent of this choice of the order parameter, and the discussion 
here is intended mainly to avoid any unnecessary vagueness that tends to lead to concerns, 
such as the very existence of an appropriate order parameter which is finite on both black- 
hole and black-string. 

For concreteness, we take the compactification manifold to be X = throughout this 
subsection. Far away from the black object the leading behavior of the radial coordinate r 
is well-defined (by comparison with the flat geometry) and thus, as usual, it is possible to 
read the (ADM) mass of the object by the asymptotic r behavior of the metric functions. 
One such asymptotic constant can be measured from the fall-off of gu and for a spherical 
hole in a flat (and topologically trivial) background this would be the only independent 
asymptotic constant, and it would be proportional to the mass. Here there is one more 
asymptotic constant: the metric becomes z-independent (z-dependent modes are massive 
from the lower dimensional point of view and hence they decay like exp(— 27rn r/L), n £ Wj) 
and thus it is sensible to perform a dimensional reduction asymptotically. After dimensional 
reduction g^z, the size of the extra dimension, turns into a scalar field. Thus we expect 
two asymptotic charges - the mass and the scalar charge. (The latter is non-conserved, but 
is conventionally called "a charge", presumably since it is the coefficient for the leading 



^Except for subsection |4.3| where a simulated time evolution is described, and subsections 2.3 6.2 where 
the open questions are discussed. 

^"We use the notation A for the general discussion of an order parameter, to distinguish it from the closely 
related perturbation parameter around uniform strings, which we denote by A, that will be introduced later 
and which also satisfies that A = if and only if the string is uniform. 
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l/r*^"^ asymptotic fall-off of a field just like the electric charge can be read from the fall-off 
of the electro-static potential. It is precisely this property of being the leading term in the 
asymptotic region which we are interested in.) 

One can define the asymptotic charges from either the higher dimensional or from 
the dimensionally reduced points of view. In the higher dimension the metric defines two 
asymptotic constants 



11 



9 



5^ = 1 - 
= 1 + 



2a 
2b 

rd—'i 



(3.1) 



From the lower dimension point of view the definition for b conforms with the standard 
definition of a scalar charge: One defines the scalar field from the Qzz component of the 
metric through e^* := Qzz, and the scalar charge A$ through the asymptotic behavior of 
$ as $ = K^/r'^~^ (we are not careful here to fix constants in any particular way). Thus 
from (3J) we see that A$ = b. 

We identify the total mass from the dimensionally reduced metric which is gotten by 
a Weyl rescaling of the metric gu — > gtt = Qzz^^ 9tt and therefore a = a — 6/(D — 3), 



where a is defined by 
get 



-gtt 



1 



G,M 



2a/r . Identifying 2 a with 



d-3 



and using 



ivr 



'-[{D-3)a-b] . 



we 



(3.2) 



In addition to the mass, the asymptotic constants a, b can be used to express another 
physical charge, the tension r. In a satisfying analogy with the well-known first law of gas 
thermodynamics 

dE = TdS-PdV, (3.3) 

where E, T, S, P, V are the energy, temperature, entropy, pressure and volume, respec- 
tively, r is defined here through 



dM = TdS + TdL . 



(3.4) 



Namely, the tension is defined to be the thermodynamic conjugate to L, the size of the 
extra dimension (see ^] for earlier and equivalent definitions of tension). 

The thermodynamic charges are related to the asymptotic constants a, b through 



M 


_ ^D~3 


tL 


Svr 



D-3 -1 
1 -{D-3) 





a 




b 



(3.5) 



This relation may derived either through the thermodynamic definitions or from the "method 
of equivalent sources" as we proceed to explain. The thermodynamic definition is fully 
specified by the gravitational (Gibbons-Hawking) action / = —PF where /3 is the inverse 



Exactly two independent ones as discussed above. 
^^This is our own notation for this known method 



see for instance 



but we shall not attempt 



complete referencing for it. 
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Uniform string 


Small black 


hole 


Scalar charge 





GdM/{D - 


-3) 


Tension 


ML-^/iD - 3) 








Table 2: Values for some possible order parameters for both "extreme" phases 



temperature and F is the free energy. Alternatively, in the "method of equivalent sources" 
one imagines that the asymptotic fields were generated by a weak stress-energy source and 
uses the linearized equations to infer the integrated stress-energy charges from the metric 
asymptotics. Note that the current expression for the mass ( |3.5| ) coincides with the mass 
read off the dimensionally reduced metric (p.2|). 

Let us gain some insight into the behavior of the tension and the scalar charge. For the 
uniform string metric (2!3|,2.4) 6 = from its definition (3J), and hence t L = M/{D — 3) 
from (|3.5|). For the small black hole, on the other hand, one finds from the Newtonian 



approximation ( 4.30 , 4.16 ) that r = and b = a/{D — 3). More precisely, to leading order 
the tension is proportional to M^: tL = {D - 3)Q{D - 3)/2 {po/L)^-^M |3^. See 
table ^ for a summary. 

Inverting (^) we obtain 



a 


Svr 


D-3 -1 




M 


b 


Qd~3{D-4) {D-2) 


1 -{D-3) 




tL 



(3.6) 



Looking at the expression for the scalar charge b we see that the mass tends to increase it, 
namely mass "wants to generate more space" for itself, while tension "wants" to contract 
the extra dimension. Thus we may say that for the uniform string (b = 0) the tension has 
exactly the correct value to cancel the tendency of mass to expand the extra dimension. 

In fact we empirically find that for all black hole and black string solutions their 
6, T parameters lie between the uniform string and the small BH. This is only partly 
understood. Positive tension r > 0, was proven in analogy with the positive mass theorem 
1 49, However, it is not clear so far why b > holds. Actually, when one considers 
also bubbles (for instance ||5l|) then b is no longer positive. While b > would correspond 
to the bound t L/M < 1/{D — 3) it was argued in |4^ (see also |52|) that the completely 
general bound is higher, namely t L/M < (D — 3). This bound is set by the bubble and is 
consistent with the Strong Energy Condition Tqq — 1/(L) — 2) qqq > 0. 

We may now define the order parameter. Since b is zero exactly for the uniform string 
we can use some multiple of it. The natural choice is a dimensionless scalar charge, being 
either b/{GdM) for the micro-canonical ensemble or b/f3'^~^ for the canonical ensemble. 
The dimensionless scalar charge has the additional advantage of placing all phases at finite 
values: not only is the uniform string at 6 = but also the small black hole is at finite 
value, namely b/Ai = 1/(D — 3), as can be seen from table ^. Alternatively, one may choose 



The gravitational action is defined in 3.13 



and discussed around it. 



The metric signature convention is "mostly plus" 
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a dimensionless tension, such as t/M, as an order parameter that vanishes not for uniform 
strings, but rather for black holes, and is closely related to h and M ( |3.5[ ). 

Here we note that the differential form of the first law of black hole thermodynamics 
( |3.4| ) may be integrated using the scale invariance of GR. Namely, when one scales the 
lengths L — > e° L the mass and area (entropy) scale according to M ^ g(D-3)a g _^ 
g(_D-2)a Taking the differentials for these transformation with respect to a at a = 1 and 



substituting into (3^) yields 

{D-3)M = {D-2)TS + tL , (3.7) 

which is a useful formula known as "the integrated first law" or "Smarr's formula" (shown 
in the current context in 1 46 , ^ ) . 

An interesting property of a phase diagram with this order parameter is that intersec- 



tions in the phase diagram are constrained due to the first law [46|. 



3.2 Order of phase transition 

In general, one of the basic properties of any phase transition is its order. In this subsection 
we first review black hole thermodynamics in the Gibbons-Hawking formalism and the 
general Landau-Ginzburg theory and then we summarize the results for the system under 
study. 

The gravitational free energy. We will use the standard semi-classical^^ Gibbons- 
Hawking gravitational free energy given by the gravitational action I{g^v'-, P) = —PF 
(to be described below) evaluated on a "Euclidean section" of the metric. Since our solu- 
tions are static, it is straightforward to obtain a "Euclidean section" simply by taking 
the transformation t ^ it. In a standard way, requiring the absence of conical singularities 
at the horizon fixes the period of Euclidean time to be /? = Itt/k, where k is the surface 
gravity which is constant over the horizon by the zeroth law.^'' 

The gravitational action / is given by the standard Einstein-Hilbert action with an 
additional boundary term Iq 

I = Ieh + Id (3.8) 

such that / is stationary on solutions with respect to variations of the metric which preserve 
the boundary metric |53, 54 1. 



^^From a practical point of view all the computations with this action are classical, h enters only in 
the dictionary between the variables of F such as (5, A (the periodicity of Euclidean time and the area), 
and the thermodynamic variables such as T~^,S (the inverse temperature and the entropy). For example 
S = A/(4Gft). 

"Static" means by convention "non-rotating and time- independent" or more formally, invariance not 
only under time translations but also under time reversal, so that gu = 0, i t. Equivalently, there exist 
hypersurfaces such that the Killing vector field dt is orthogonal to them, namely the t — const hyp ersurf aces. 

^^The zeroth law is derived by imposing the constraint Gni ~ 0, where G is the Einstein tensor, n is the 
coordinate normal to the horizon and axe the coordinates tangent to the horizon, excluding t, the time. 
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where R is the Ricci scalar and Iq can be defined by either of the fohowing equivalent 
definitions^^ 



SirGIa := / K (3.10) 
JdM 

where K is the trace of the second fundamental form for the embedding of the bound- 
ary in the manifold. 



8TrGl9:=dnV9, (3.11) 

the derivative of the boundary {D — l)-volume with respect to a (proper length) shift 
in the normal direction. 

• Iq is the boundary term obtained through integration by parts of Ieh such that 
/ contains only first derivatives of the metric and no second derivatives, namely 
{dg){dg) G /, gd^g ^ / where g denotes here a generic metric element. 

In asymptotically flat spaces Iq as defined above (|3.10| ) diverges and must be regular- 



ized. The standard regularization |53] is done by measuring I relative to flat space. One 



chooses a large cutoff r = R, where in our case the boundary is S^~^ x S]^^^^ x S^^^^, 
where L{R), (5{R) are the periods of the 2, t directions, respectively"*^^. Next one subtract 
the action of a flat space with the same boundary, which in our case is R*^"^ x S]^^^^ x S^j-^-j, 
namely 

SttG/o = dr (nd-2 r^-^ L{R) (3{R)) \^ = {d- 2) f^d_2 L{R) P{R)R'^-^ . (3.12) 

Finally one takes R ^ 00. 

Combining this with (3.S, 3.9 , |3.1Cl|) , the free energy F = F{g^y](3) is finally given by 



^ = ^=T^f ^+^/ (^-^0) (3.13) 



which includes a bulk integral over R, the Ricci scalar and a boundary integral over K — Kq, 
where K is the trace of the second fundamental form on the boundary, and Kq is the same 
quantity for the reference flat space geometry. 

Landau-Ginzburg theory. In a phase transition some derivative of the free energy is 
discontinuous and goes through a jump. The order of a phase transition is defined to be the 
order of this derivative. A first order transition is between two phases which are separated in 
configuration space and hence have different entropies (and other thermodynamic variables) 
and are therefore exothermic (involving latent heat) while for second order and higher the 
phases are continuously connected and there is no finite release of energy. 



|53[ uses definition ( p.lCj ), and the equivalence with definition (3.11) is used in a computation. I beheve 
that the third and last definition is implied by the relation between the boundary conditions and the action. 
^^For large r the black object metric is virtually z independent, as we discuss above 3.1. 
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The Landau-Ginzburg theory of phase transitions |55] tells us how to infer the order of 
the transition from the local behavior of the free energy near the critical point F = F(X; //), 
where the A variables parameterize (a selected subset of) the configuration space, and 
may act as order parameters, and the fi variables are the control parameters, for instance 
the temperature. In our case, the dimensionless control parameter may be chosen as 
fi = fi/3 (X (3/ L. Geometrically it is the ratio of the two asymptotic periods, while physically 
it is a dimensionless inverse temperature. Since the control parameter is essentially the 
temperature, we interpret the thermodynamics as taking place in the canonical ensemble, 
and accordingly, the name "free energy" is fitting. The relevant configuration variable is A, 
the amplitude of the marginally tachyonic GL mode ( 2.11| ), so symbolically the perturbation 
is 

~ A exp(i k z) hcPY ■ (3-14) 

A generic phase transition is of first order as depicted in figure ^. ^'^ However, in 
our case F possesses a certain parity symmetry (which is non-generic) which opens the 



possibility for higher order transitions, as we proceed to explain. We note from (3.14) that 
A is complex and its phase is related to translations in the z direction. Since the action is 
invariant under z-translations we have 

F = F{^l,\X\^) , (3.15) 

and the non-uniform phase spontaneously breaks this symmetry. From now on, without 
loss of generality, we consider A to be real and omit the absolute value notation, so that 
F is an even function, as claimed. This corresponds to fixing the z-translations in ( 3.14| ) 



through /i ~ A cos(A; z) hcpY with a real A, and the sign reversal now amounts to shifting 
z by half a period. 

Note that A can be related to the order parameter defined in the previous section, the 
scalar charge b. The latter being invariant under z-translations must be a function of |Ap 
as well. Since they both vanish for the uniform string we conclude that 

6oc|A|2 (3.16) 

(there is a genericity assumption made here which is confirmed by calculations). 

Having a marginally tachyonic mode appear at some critical value fic means that the 
quadratic term in A has a zero at fic, namely 

F = Fi{i2-fic)X^ + ... (3.17) 

where Fi is some positive constant (assuming the phase is stable for fi > Hc as is our case). 
Now the Landau-Ginzburg theory suggests to expand the free energy at fic to higher orders 
in A and to test whether the free energy has a minimum or not^^. If F has a minimum at 



^°See a discussion of the Hawking-Page transition in subsection |3.3| . 

^'^In the general case, where F may also have negative modes, the precise criterion for whether the 
transition is first order is whether the fi of the new phase increases. Note that in the expansion one may 
need to incorporate the back-reaction. These issues are further discussed below 



-18- 




Figure 7: The most generic phase transition - first order resulting from a non-zero cubic term in 
F. On the right we see a sequence of free energy functions, parameterized by /i, with their extrema 
(phases) highlighted. On the left the phases are extracted into a phase diagram. Stable (unstable) 
phases are denoted by solid (dashed) lines. 





(a) 



(b) 



Figure 8: For an even free energy the order of transition could be first or higher depending on a 
sign. This figure shows a second order transition (or higher) since the free energy has a minimum 
at the critical point. Same conventions as in figure ^. 




Figure 9: An even free energy results in a first order transition if it has a maximum (or more 
precisely a non-minimum) at the critical point. Same conventions as in figure ^ 



He as in figure then for ^ < jic two stable minima are created close by at small A, the 
system will continuously evolve into these new phases and the transition must be second 



-19- 



order (or higher). If however F has a "direction of descent" at /ic, as in figure that 
means that since F must be bounded from below there must be some other minimum at 
some finite value of A, whose free energy is lower than the A = phase. Therefore the 
system underwent already a first order transition at some higher value of where the free 
energies of both phases were equal. 

Therefore Landau-Ginzburg theory instructs us (in our case, where F is even) to 
expand 

F = Fi (/i - /xe) A^ + F2 A^ + . . . (3.18) 

and to determine the sign of F2, the quartic coefficient at /ic- Positive (negative) F2 implies 
a second (first) order transition. Of course if F2 happens to vanish one needs to compute 
higher orders, but this did not happen in this system. 

Actually, so far we neglected all other modes except for the GL mode. When these 
are brought into consideration one finds it is required to compute first the (quadratic) 
back-reaction of the GL mode and incorporate it into the computation of the quartic term 
in the action. The reason for taking the back-reaction is that what we really want to know 
is whether the emergent phase from the critical point goes up or down in ^ in the phase 
diagram, indicating a first or second order transition. To see that we denote the extra 
modes by and expand the free energy as follows 

F(A, y';^i) = Fi{fi- /x,) A^ + F2X A^ + F2i X^y' + K.j y' y^ + . . . (3.19) 

The equation of motion for the is = = '^Kijy^ + F2i A^ = 0. We denote the 

solution by = A^ 5*, where -B* stands for back-reaction. The equation of motion for A is 
= ^F = 2Fi(/i — ^c)A + 4 F2\ A'^ -I- 2 F2i Xy\ After diving by A and using the equations 
of motion for we find that -2Fi(^ - fi^) = 4:X^ (-F2A - Kij B^) . Namely, in order to 
determine sgn(/x — ^c) and from that the order of the transition we need to compute 

F2 := F2X - Kij B' B^ (3.20) 

where B^ is the back-reaction that solves the equations of motion. Note that F2 is 
computed without deviating from criticality fi = fic- 

Results. The determination of the order was first carried out by Gubser |^] for the back- 
ground M^'^ X (see also Wiseman's improvement (5^). Part of the original motivation 
there was to find a second order phase transition and hence together with it a branch of 
stable non-uniform strings emanating from the GL point, such as those predicted by p3| . 
However, the transition was found to be first order. Sorkin generalized the method to 
the backgrounds M^^^'^ x and found a surprising critical dimension 

^^A non-standard term which we use to mean "some direction in configuration space where F decreases". 

^^The index i runs over all the extra modes, namely the additional configuration variables. In our case 
it should really be the continuous variable r and sums should be replaced by integrals, but we keep this 
notation for conceptual clarity. We expended F up to quartic order in A, incorporating the fact that will 
be second order in A. 
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Figure 10: The first correction to the mass rji and to entropy (relative to a uniform string with 
the same mass) (72 for the non-uniform string branch emanating from the GL point for backgrounds 
■^D-2,1 ^ gi _ gj j-Qj. pj-gcigg definitions. From our general discussion we know that there is 
only one sign to be determined (that of F2), and indeed the signs of rji, (T2 are correlated. At the 
critical dimension D* = 13.5 there is a change of sign indicating the transition becomes second 
order for higher D. Adapted with permission from B3]. 



L> < 13 


first order 


D>U 


second order . 







The critical dimension D* = "13.5" is demonstrated in figure Kudoh and Miyamoto 
1 58] observed that the critical dimension depends on the ensemble: Sor kin's critical dimen- 
sion holds for the micro-canonical ensemble, while in the canonical ensemble the critical 



dimension is lower by one. In subsection 6.1, "results", we comment on the relation of 
these results with the predictions of pSj . 

Due to the importance of the critical dimension it is a good idea to develop some 
intuition about it. First, as discussed around equation ( p. 121) for high D the critical GL 
string is quite "fat" and we expect that the string will not decay (directly) into the black 
hole, which would be "too big to fit" inside the extra dimension (see also for a similar 
argument involving ^5, which is defined below). Another indicator comes from comparing 
UGL with us the value of the approximate equal area point: SstilJ'S) = S'sh (/^s) where 
Sst is the entropy of the uniform string and Sbh is BH area approximated by the small 
black hole expression. For a first order transition we expect /X5 > ficL but that holds only 
for D < 12.5. 



Of course dimensions are integral, and the notation means only that the change in the order happens 
between D = 13 and D = 14. 

■^^A black hole embedded in a heat bath. 
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Method. Gubser's method [^] was to perturbatively follow the branch of non-uniform 
strings emanating from the GL critical point. The perturbation parameter is A, the ampli- 
tude for the GL mode ( |3.14| ), and in order to determine the order of the transition is was 
necessary to compute some metric functions up to the third order in the perturbation pa- 
rameter. As we mentioned above, to determine the order there is a somewhat more efficient 
method, namely it is actually enough to compute only the second order back-reaction and 
substitute into the quartic part of the action (see appendix A in [^] and |]l8[). However, 
the longer computation naturally yields additional results not included in the shorter one. 

In practice, when computing the back-reaction we consider a continuum of modes, 
and therefore the discrete index i in ( |3.19| ) and the discussion around it is replaced by 
the continuous variable r, and the "field"?/ is replaced by the all the fields in the problem. 
Moreover, inverting the linear operator K to solve for y means obtaining the fields by 
solving a second order ODE (obtained from linearizing the equations of motion) with a 
source term quadratic in the perturbation (the first order mode). 

3.3 Morse theory 

When one turns to consider the question of the end-point for decay, or more generally 
of finding the phase diagram of all static phases, one is at first discouraged by the lack 
of any knowledge regarding the non-uniform strings and the big black holes. The most 
interesting question is to determine the qualitative features of the phase diagram. For that 
purpose one needs qualitative tools, and one such tool was given in |4^ ] under the name 
"Morse theory", fulfilling the intuition that generically a phase persists as the system is 
"deformed" by changing a parameter, and that the disappearance of a phase should require 
some special circumstances that are worth elucidating. 

Indeed solutions of the Einstein equations are extrema of the gravitational action in the 
space of metrics, and as such are generically stable under perturbations. The topological 
theory of extrema of functions is well-known and is called "Morse theory" and it includes 
the specification of the allowed transitions. 

For other qualitative tools in the study of thermodynamics in the astrophysical context 
of self-gravitating systems see the review |6^, and especially the closely related Poincare 
method to determine the perturbative stability of phases just by looking at a certain kind 
of a phase diagram. 

Actually, there is a subtlety in the identification of extrema of the action with solutions 
of the equations of motion due to gauge (diffeomorphism) redundancy, which we would like 
to mention. It is certainly true that solutions are extrema of the action. However we wish to 
consider the action as a function of metrics up to gauge invariance, namely as a function of 
the gauge-fixed metric. Therefore the extremum equations should be supplemented by the 
gauge-fixing constraints. Nevertheless, in this case it was found that the constraints are 
actually implied by the extremum requirement through a combination of properly chosen 

^''Some readers may be familiar with the way Morse theory measures global properties of manifolds 
(Homology), but here we need a different aspect of the theory - local invariants of extrema under deformation 
of the function. 
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boundary conditions and the constraints' Bianchi identities. It could be that this is true 
more generahy. 



A lightning review of Morse theory (see section (3.2) of [^] for a somewhat longer 
introduction). For functions of one variable the ways in which an extrema can disappear 
are clear 

• Annihilation. A maximum and a minimum can coalesce under continuous deforma- 
tion and disappear into a monotonous function. We call this the basic Id vertex 
of "annihilation" - see figure 0. 

• Run-away. An extrema of a function f{x) may run away to infinity either in x or 
in / during a finite range of deformation parameter. In this paper we shall find 
"annihilation" explanations for changing phases, and thus we will not need to resort 
to "run-away" explanations, even though they are certainly a logical possibility. 

When one considers a function of several variables, / = /(x), one may get a simple 
generalization of the basic Id vertex by adding n spectating negative directions as well as 
N-n-l positive directions, namely /(f) = /o(xo; /u) - Y17=i h + Yl!j=n^\ fj ^l-. where 
/o(xo; [i) is a Id function such as the one depicted in figure ^, and /j, \ < i < N — \ are 
positive constants. Now, the minimum in figure |^ turns into an extremum with n negative 
modes, while the maximum has n -|- 1 negative directions. Therefore we obtain, what we 
call "the basic vertex" 



Basic vertex: two extrema with n and n + 1 negative directions may annihilate. 



While for generic extrema the Hessian is non-degenerate and n is well-defined, one may 
wish to know the rules for more general vertices. Indeed, Morse theory can be phrased 
as saying that the most general vertex is a coincidence of several of these basic vertices, 
thereby justifying our use of the adjective "basic". 

The conclusion and some reservations. We see that a stable phase in = 0) is allowed 
to disappear at the expense of "annihilating" with an unstable phase with one negative 
mode (n = 1). The latter in turn can disappear by annihilating against either n = or 



n = 2 phases and so one. This is exactly the ''''phase conservation rule" |44] which we were 
seeking. It sets a strong qualitative constraint on the existence of phases. However, the 
price to be paid is that all phases must now be mapped out, not only the stable ones. 

As shown in figures this rule determines the stability of the non-uniform string 
emanating from the GL point. One may ask where this phase might end. From the phase 
conservation rule we conclude that the simplest way to satisfy it, without requiring any 
additional phases, is that the non-uniform string phase would annihilate against the black 



^'^Here we introduce the term "vertex" to mean the event when two or more extrema of a function coincide 
as a deformation parameter is varied. Such an event looks hke a collision of phases in a phase diagram and 
the name comes from the analogy with the Feynman diagram vertex at the collision point between two or 
more particle world-lines. 
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hole phase |44] at a point on the phase diagram that we call "the merger". Taking into 
account the existence of a critical dimension we find that depending on the dimension we 
get two possible behaviors: for 5 < D < D* the non-uniform black string, which is unstable 
for small non-uniformity, annihilates with the black hole (which is stable when it is small) 
while for D > D* a stable non-uniform string transforms into a stable BH, and D* = 12.5 
(13.5) in the canonical (micro-canonical ensemble). This is our main conclusion from Morse 
theory and we stress again that it seems to be the simplest scenario, but others cannot be 
excluded. 

The rigor of the prediction of a phase merger, even if intuitive and clear, is questionable 
due to the following observation. In the next section we will see that the Euclidean versions 
of the black hole and the string have different topologies and hence their metrics would be 
expected to live in different, disconnected spaces of metrics, and it wouldn't make sense 
for phases to move from one space to the other. Nevertheless, we shall take the prediction 
above seriously and look whether these two spaces of metrics are in some sense glued 
together. Indeed we shall find a continuous transition (and in finite distance) between the 
two spaces, which we view as an important confirmation for the consistency of the picture. 
However, the way in which the spaces are glued is still poorly understood, and the gluing 
may very well be non-smooth as well as involve the infinite dimensionality of the space of 
metrics in an essential way, in which case the validity of the Morse theory argument is not 
self-evident. At this point, I consider the conclusion above to be essentially correct and 
justified if not rigorously a priori then a posteriori by the agreement of the predicted and 
the numerically computed phase diagrams. 

An example: the Hawking-Page transition. The reader familiar with the phase 
transition between thermal Anti-de-Sitter (AdS) and large black holes in AdS, known as 



the "Hawking-Page transition" |61], may benefit from applying the Morse theory ideas 
above to that context. 

In the Hawking-Page transition in its canonical ensemble setting, one considers a space- 
time with a negative cosmological constant A whose boundary is that same as that of 
thermal AdS, namely a large spatial sphere of Radius R times a compact circle of Euclidean 
time of size P\/'^^R. The dimensionless control parameter is /3a := P V— A, and as usual 
(3 /his the inverse temperature. In 4d Hawking and Page found that several phases exist 
in this system as follows. For Pk > (3o = the only phase which exists is thermal AdS, 
AdS filled with a thermal gas of radiation. For do > dh > Pi = —k- two additional phases 
show up, the small and large AdS black holes (inside a thermal bath), the large (small) 
black hole being thermodynamically stable (unstable) due to its positive (negative) specific 
heat. However, the black holes' free energy is inferior to that of thermal AdS. Finally at 
/3a < /3i the large black hole dominates. 

In the current context, this phase transition, which is a first order phase transition 
is described in figure |ll| (compare with the general first order transition of figure 0) . For 
P > there is a single minimum for the free energy which is thermal AdS. At /? = /?o two 



Another transition is expected at /32 oc h}'''^ which is of quantum nature and will not be discussed here. 
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Figure 11: The free energy and phase diagram for the Hawkmg-Page transition in AdS space, 
which is a first order phase transition. At (3 = (3^ the unstable smaU black hole phase meets the 
stable large black hole phase and they "annihilate" according to our "basic vertex" rule. 



phases are "pair created" (or "annihilated" if one approached Pq from below) through our 
"basic vertex" : a stable large black hole and an unstable small black hole, but thermal AdS 
still has lower free energy. Here we are assuming the Gubser-Mitra Correlated Stability 
Conjecture pi, 37] to infer the perturbative stability or instability from the sign of the 



specific heat and the associated thermodynamic stability of the black holes. Then at 
P = Pi a first order phase transition occurs when the free energies of thermal AdS and the 
large BH become equal. Finally, for smaller P the large BH dominates. 

3.4 Merger point 

In the last subsection we saw how Morse theory makes it plausible that the non-uniform 
string phase merges with the BH phase. We shall first encounter a problem for this picture, 
namely topological differences between the two phases. However, we are familiar with some 
topology changes such as the flop and the conifold, and one of the surprising results of the 
research on this system is the emergence of a novel type of topology change, called the 



"merger" transition in |44|. 



A topology change. Intuitively the transition from black string to black hole involves a 
region where the horizon becomes thinner and thinner as a parameter is changed until it 
pinches and the horizon topology changes. We call this region "the waist" and this process 
is described in the upper row of figure 12 using the (r, z) coordinates defined in figure |2|. 



It is important to remember that all metrics under consideration are static and that they 
change as we change an external parameter, not time. Since the metrics are static we may 
as well consider their Euclidean versions (this point was discussed in the second paragraph 
of section IsT 
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We shall now demonstrate that this merger transition involves a local topology change 
of the Euclidean manifold. Let us zoom in around a very thin waist, whereby all the scales 
of the problem such as GM and L are very large and irrelevant, identifying what may be 
called the local geometry at the waist. Consider a co-dimension 1 surface within the local 
geometry but far away from the waist, such as the one denoted by a dashed line in all three 
geometries in the top row of figure Actually, this surface is the asymptotic boundary 

of the local geometry. 

The topology of this asymptotic boundary is given by S^~'^ x ^ ^ as we proceed to 
explain. The angular piece, S^~^, is obvious while the S^^^^ requires explanation. One 
should remember that in our figures, such as figure we suppress not only the angular 
coordinates but also the time variable. In the Euclidean continuation the Euclidean time 
must be periodic in order to avoid a conical singularity at the horizon, and moreover, the 
proper size of this circle vanishes at the horizon [gu = at the horizon). Thus the circle 
fibration of Euclidean time over an interval (the dashed line in figure [l2| ) , such that the 
fibred circle shrinks on the edges exactly produces a topological S^. This is just like the 
fibration of the surface of the earth by latitude lines, which shrink at the poles (see figure 
Vc 



Naturally, the topology of this boundary surface is constant as local changes occur 
near the waist. Next we notice that in the black-hole phase the S^~^ is contractible (onto 
the exposed r = axis), while in the string phase the S^^t is contractible. Therefore the 
local topology of (Euclidean) spacetime is changing, not only the horizon topology. Thus, 
the topology change can be modelled by a "pyramid" (just like the conifold transition) - see 



the lower row of figure 12: the rectangular basis of the pyramid denotes the asymptotic 



boundary S^~^ x ^ ^ where each edge represents one of the sphere factors, and the 
pyramid's truncated apex encodes which one of the spheres is contractible and which one 
remains of finite size and is non-contractible. By the nature of topology, in order to change 
it there must be at least one singular solution along the way (with at least one singular 
point). The simplest possibility, (which is also realized in the conifold) is to assume that 
the the singular geometry is the cone over S^^^ x S^. 

It is suggestive that the local space-time topology change is also accompanied by a 
change in the global topology. I believe this is the case for the following reason. The black 
string solution is contractible to S^~^ x (after contracting the r,t cigar), and thus its 
elementary non-contractible cycles are S^~^ and S^. For the black hole on the other hand 
the non-contractible cycles seem to be S^^"^, S^ j and S^: the S^g~^ is the horizon, the 
is the compact z direction and the ^ is the r = axis connecting the poles together 
with the time fibration. It is seen that there are several topological differences between the 
geometries, for instance, the black hole has a 2d topologically non-trivial cycle while the 
string does not. 



^^In other words, if we denote by pw the radius of the (minimal) angular sphere {S^~^) at the waist, 
then we wish to consider the fixed p surface where p,„ <^ p <^ L. 
^"which is topologically trivial since ga — for all i ^ t. 
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More formally, contracts to a point. 
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Figure 12: The merger transition. Shaded regions are inside the horizon and the dashed hne is a 
boundary far away. The singular configuration is a cone over x S^. 




Figure 13: An iUustration of the fibration of an S^: circles fibred over an interval such that the 
fiber shrinks at the edges (left) are equivalent to a surface of a sphere (right). In this way the 
dashed intervals in the top row of figure 12 actually represent spheres, once the circles of Euclidean 
time are taken into account. 



Cones. Figure 12 encodes the topological nature of the merger. But is it really true that 
this can be realized with Ricci-flat metrics? 

It is easy to write down a Ricci flat metric for the singular solution, the cone (see the 
middle of figure |l^ . Actually this can be done for a somewhat more general cone, the cone 
over S™ X S" with no additional "cost". The metric is 

ds^ = dp^ + -jj-^ [{m - 1) d^l^ + (n - 1) , (3.21) 

where the p coordinate measures the distance from the tip of the cone, D = m + n + 1 
is as usual the total spacetime dimension and the constant pre-factors are essential for 
Ricci-flatness. Note that p = is the singular tip of the cone, unless m = (or n = 0) 
when it becomes the smooth origin of M'^ in spherical coordinates. 

In order to exhibit "smooth cone" metrics which approach the singular cone from both 
sides of the transition (see right and left portions of figure 12) one may use the following 
ansatz 

ds"^ = dp'^ + dnlm + ''(^) dnl„ (3.22) 

with boundary conditions at p ^ 

a{p = 0) = ao 

b{p) = log(p) (3.23) 
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such that S™' becomes non-contractible while S" joins with p to make a smooth neighbor- 
hood of the origin of M"^^. 

Using ( [A.ip for the Ricci scalar of a fibration and after integration by parts one gets 
the action for a, h 



S= / dpe""'^^'' L(m-l)(a'2 + e-2'^)+n(n-l)(6'2 + e-2*) + 2mna'6' 

Einstein's equations are 



(3.24) 



// / 1 \ -2 a / 2 / T / 

a =[171—1)6 —ma —nab 

I," I -i\ —2b t/ 2 / T/ 

=(n — Ije — no — ma b 

= m{a' ^ + a") + n{b' ^ + b"). (3.25) 

These equations are very similar to the equations encountered in the Belinskii-Khalatnikov- 
Lifshitz (BKL) analysis of the approach to a space-like singularity (see the recent excellent 
"Cosmological Billiard" review Q). Although the general (and often singular) qualitative 
behavior of these equations as p — > for arbitrary initial conditions was not obtained in 
|44],^^ it was checked that for the boundary condition ( 3.23| ) a{p), b{p) — log{p) can be 



expanded in a Taylor series around p = and the recurrence equations for the Taylor 
coefficients could be solved without encountering an obstruction. 

Once a single smoothed cone solution is available, constructing a full family that ap- 



proaches the singular cone is a matter of simply rescaling it. As illustrated in figure |14], since 
away from the smoothed tip the smoothed cone asymptotes to a cone, a geometry which is 
scale invariant, then after rescaling there is a natural way to identify the asymptotic cones, 
thereby specifying the way to take the limit over the family of rescaled metrics. 

Altogether we succeeded in realizing the local topology change encoded in figure |l^ by 
a family of Ricci-flat metrics. I consider this non-trivial property to be strong evidence for 
the merger picture as presented in subsection \i.3\ 



I would like to stress some of the assumptions involved in locally modelling by cones 

• The singular solution has a single singular point. 

• The singular solution is continuously self similar (CSS). 

Both assumptions are reasonable and minimal: there could be more than a single sin- 
gularity, but there is at least one, and the local singular solution must forget the "long 
distance" scales and hence it would be scale-free, and the simplest way to obtain that is 
if the solution is self-similar. Continuous self-similarity would be the simplest possibility 
and in |44| it was found to lead to a local isometry enhancement where the x coordinate 

■^^One can prove useful theorems for the evolution of the volume factor U = 6"""+"* using the geometric 



analogue of the "c-theorem" (I thank J. Maldacena for pointing this). From (3.25) we have (log(l7))" — 
— m a'^ — nb'^ < Hence if U (or equivalently log(l7)) is somewhere decreasing it must continue to decrease 
monotonically. At the same time U" — (m(m — 1) e"^" + n(n — 1) e~^'') U > which together with the 
previous result guarantees that U is monotonic. 
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Figure 14: A single smoothed cone solution can be scaled down to provide a continuous family of 
metrics which approach the cone. 



(see figure conspires with t to make a round S^. Indeed, a Numeric study [S2|, limited 
by numerical resolution, found consistent evidence for this cone structure. More recently 
it was claimed in |6^] that in a certain range of dimensions (5 < D < 10) continuous self- 
similarity (CSS) is in fact (spontaneously) broken into discrete self-similarity (DSS) which 
requires to replace the cones by certain "wiggly cones" as local models of the merger point. 

Tachyons on cones. It turns out that the cones may have tachyons and that their 
existence surprisingly depends on a critical dimension -D^g^ggj. = 10 as we proceed to show 



- see also |44| and references therein including |64|. 

The dangerous mode is a function e(p) which inflates slightly one of the spheres while 
shrinking the other. The ansatz for the perturbation is 



D 



m 



1) ^Z™ + (n - 1) ^/'^ 



(3.26) 



A priori one could start with two separate scale functions, one for each sphere, but the 
Gpp = constraint relates them as above (after ignoring the trivial perturbation which 
represents /o-translation) ~ for more details see |Q eq. (6.6) and below. 
The quadratic part of the action, disregarding an overall constant is 



dp 



(3.27) 



and through the change of variables p ^ = [D — 2) p^ ^ it can be recast to have a canonical 
kinetic term 

2 



dp 



-e' 2 



The equation of motion for e, namely the zero mode equation, is 



(3.28) 



+ V{p) 
V{p) 







(Z)-2)p2 



(3.29) 
(3.30) 
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The solutions are 



^1,2 = ^^ . (3.31) 



The expression (3.31) for si^2; the characteristic exponents, reveals a critical dimension 

m 

D* = 10 (3.32) 

such that for D < D* s are complex while for D > D* they are real. Complex charac- 
teristic exponents (for D < D*) imply that (the real part of) the zero mode has infinitely 
many nodes (zeroes) equally spaced in log(p) with "log-period" Ag = 27r/9(s). The pres- 
ence of infinitely many nodes for the zero mode implies the presence of infinitely many 
tachyons, just like the number of nodes of the zero-energy solution to a Schrodinger equa- 
tion counts the number of negative energy states. In [^] these log-periodicity and tachyons 
(for D < D*) were interpreted as an indication for a spontaneous breaking of continuous 
self-similarity (CSS) into discrete self-similarity (DSS). For D > D* on the other hand, 
9(s) = and there are no nodes nor tachyons. 

Note that D* is the total spacetime dimension, and it is independent of m and n 
separately. When we wish to distinguish this critical dimension from others we shall denote 
it by -C^merger- 

One can view the zero-mode equation ( |3.29| ) in a wider context by considering the 
eigenvalue problem 

[-9^ + V{p)]e = Xe (3.33) 

where V is the same as in eq. ( |3.30| ), A is the eigenvalue,^^ and we note that ( |3.33| ) reduces 
to ( |3.29| ) upon setting A = 0. 

The eigenvalue problem ( |3.33| ) is in Schrodinger form, and we may apply known results. 
For potentials of the general form V = —c/r°' it is well-known that while the classical 
energy is unbounded from below, the quantum problem may have a ground state as long 
as the potential is not "too negative". For instance, for a = 1 we get the Hydrogen atom. 
More generally the spectrum is bounded from below for a < 2, while for the critical value 
a = 2, which is our concern here, the prefactor c becomes dimensionless and the potential is 
conformally invariant. Due to scale invariance the spectrum is constrained to be invariant 
under positive rescaling of the eigenvalues. Now c itself exhibits a critical value, namely 
c* = 1/4, such that the spectrum is bounded from below and actually non-negative only 
for c < c* (since the E = solution has no nodes - see for example |^5|). Equating 



2/{D — 2) = c* = 1/4 we arrive once more at D* = 10 as in (3.32). 



Some moduli space properties. The merger transition lies in finite distance in moduli 
space (actually, it would not deserve to be called a "topology change" otherwise, since it 



^^This A is unrelated to the A introduces in subsection 3.2 to represent the amplitude of the perturbation, 
and will not be used elsewhere in this review. 
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would require infinite "resources" to be reached.) The argument was unpubhshed so far 
and here it is supplemented as appendix 1^. Moreover, the appearance of a kink^^ in the 



phase diagram at merger was predicted in section (5.3) of [44|, and the numerics indeed 
seem to exhibit some sort of a kink - see the next section. 

3.5 Phase diagrams — predictions and data 

We may now assemble all the theoretical input and draw the predicted qualitative form of 



the phase diagram for various spacetime dimensions ~ see figure |I5[ These figures apply 
with small changes to both the micro-canonical ensemble, which is the usual physical setting 
where energy is conserved, and to the canonical ensemble, namely black holes immersed in 
a heat bath. 

The vertical axis is the dimensionless control parameter /i defined by either the di- 
mensionless mass ^ m ( |2.1| ) in the micro-canonical ensemble or by the dimensionless inverse 
temperature /x^ ( |2.2D in the canonical. The horizontal axis is the order parameter chosen 
in subsection to be the dimensionless scalar charge (see especially eq. (3.1) and the 
next to last paragraph in that subsection). 

The local structure around the GL point is determined by the order of the transition 
as discussed in subsection 3.2: for D < D* the GL vertex is first order |56| just like figure 



m, while for D > D* it is a second order vertex just like figure ^ D* ^ the critical dimension 
depends on the ensemble - in the micro-canonical case it is D* = "13.5" while in the 
canonical case it is D* = "12.5" p^]. 

Finally, we connect the black hole and black string phases, as suggested by the central 
conclusion of subsection |3.3| , based on Morse theory arguments and further justified by 
presenting a novel topology change in subsection 



The diagrams in figure [15| were constructed by attempting to draw the simplest diagram 
consistent with the data and assumptions in the last two paragraphs. The non-trivial nature 
of these diagrams is best illustrated by the various other possibilities that were considered 
in the literature, see for example the six scenarios in |6^, section 6. 

The critical dimension D^^^g^^ = 10 ( |3.32| ) affecting the stability of the x S" cone 
in the critical merger solution is not visible in figure According to the picture that 
emerges from |^^, there is a single branch of non-uniform strings, irrespective of D^Q^g^^- 
(and being unstable as long as the transition is first order), but for D < 10 this solution 
approaches a discretely self-similar (DSS) solution near the singularity, while for > 10 it 
approaches a continuously self-similar (CSS) solution there. While the qualitative predicted 
form is similar, the kink at merger is likely to be different. 

Another point to note is the critical point GL' in figure ^ which is there to remind us 
that each (non- uniform) solution has "harmonies" or "copies" gotten by trivially fitting 
several cycles inside L, namely replacing L ^ L/m for any m £ S+ in the solutions. In 
particular the GL point has these copies. However, these copies of the phase diagram are 
decoupled (except for their connection with the uniform string) and therefore there is no 
need to draw them. See also |6£]. 



*More formally, a discontinuity of the tangent to the phase line. 



- 31 - 



string 



GL 



•0- 



GL' \J " 




D> D' 



string 

Stable 

non-uniform 
string 



merger 




Figure 15: The predicted qualitative form of the phase diagram. The vertical axis is the control 
parameter fx and the horizontal axis is the dimensionless scalar charge, which plays the role of the 
order parameter (see the text for further detail) . Solid (blue) lines denote stable phases and dashed 
(red) lines denote unstable phases. The "+" denotes the merger point where the black- hole and 
black-string phases meet in a topology changing transition (of the Euclidean solutions). Note the 
dimensional dependence of the qualitative form due to the critical dimension. The diagram on the 
left basically appeared in Q. 



Numerical data. Clearly, the predicted phase diagrams in figure 15, were not proven here, 
but rather argued to be the simplest possibility which is consistent with certain carefully 
analyzed arguments, some of them not fully understood yet. As such it suggests to be 
tested by actually obtaining these solutions. Indeed, one of the joys of this problem is the 
feedback between theory^^ and numerics, which is in many ways like the classical feedback 
between experiment and theory which we sorely miss. 

Recently the numerical determination of the phase diagram in 6d was all but completed 
1 34] - see figure 16. A visualization of the merger transition through embedding diagrams 



is shown in figure 17. 



We see perfect agreement with the predicted diagram in figure 15 with D = 6 (which 



basically appeared already in especially regarding the prediction for a "merger" of the 



■^^The word "theory" is used here to loosely mean all the considerations that one can apply before any 
exact solutions are available 
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Figure 16: Numerical data for the phase diagram in 6d. "US'" denotes the uniform string branch, 
"NUS" denotes the non- uniform string branch (data from [^) and "BH" is the black hole (data 
from [Q). Axes conventions are the same as in figure |l^ (vertical is iip and the scalar charge is in 
units of mass rather than temperature). Adapted with permission from | |3^ . 



black hole and the string phases, and the absence of a stable non-uniform phase (predicted 
by [^Sj). We view this as a vindication of the picture presented here. Additional interesting 
features of the numeric figure are a kink at merger, perhaps related to the predicted kink 
(see the last paragraph in the previous subsection), and the location of the merger at 
roughly a local maximum on the diagram. 

These diagrams were obtained by combining data from two different simulations - one 

for previous 



for non-uniform strings |57] and one for the black holes [34|, see also 
work. More details on the challenges overcome in performing these simulations appear in 
the next section. 



4. Obtaining solutions 

In the previous section we examined the qualitative features of the phase diagram. Now 
we turn to the more quantitative aspects, those required in order to obtain solutions for 
static black objects. 

4.1 2d gravito-statics 

Counting degrees of freedom. The most general metric which is static^^ and spherically 
symmetric is 



where all functions are defined on the (r, z) plane, ds^^ is an arbitrary metric on the plane 
and since the metrics are static we might as well work with Euclidean signature. 



However, the geometry of the horizon and axis seems to merge rather smoothly. 



See also the discussion above 3.8 
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Figure 17: Embedding diagrams for a sequence of event horizons in 6d, starting with non- uniform 
strings (top, data from |^) and passing through the merger transition to black holes (bottom, 
data from X is the compact dimension and Y is the radial direction. Note that although the 

asymptotic size of the extra dimension, L, is kept fixed, the "proper" size (for embedding purposes) 
changes. Reproduced with permission from |Q. 



Altogether the problem is defined in the 2d (r, z) plane and the field content is the 
metric and two scalars A, C. That means that we can write down a 2d action for these 
fields without loosing any of the equations of motion. The action is 



167r Gd 

[R2 + {D-3){D- 4) + {D - 3){D - 4) {dCf + 2{D - 3) {dA){dC)\ (4.2) 



where R2, dV2 '■= y/g2drdz, Vtii_2 are the 2d Ricci scalar, the volume element, and the 
area of the unit sphere^*^ S"^"^ (see appendix ^ for useful formulae to determine this and 
related actions). The total number of fields is 3 + 2 = 5: 3 for the 2d metric and 2 for the 
two scalars. Two fields may be eliminated by a choice of coordinates in the plane which 
leaves us with three fields. As we proceed we shall review some of the gauges that were 
used. 



See the definition below (2 
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If we formally compute the number of "dynamic" or "physical" degrees of freedom we 
get a total of —1 + 2 = +1: —1 is for the metric degrees^^ and +2 is for the two scalars. So 
far nobody succeeded in reducing the problem to a single field, and it is not clear whether 
that is possible or not, but there is a clever ansatz due to Harmark and Obers which reduces 
the problem to two fields |6^], (see ( [4.17 , 43s| ) and thereabout) . Morally speaking, one 



may hope that the equations for the three fields, if not reducible to a single field, could 
at least be separated such that first one solves two equations for "diffeo gauge fields" and 
only then a single equation is solved for the "physical field" . 

Note that if we relax the static requirement the number of degrees of freedom increases 



see subsection 4.3 



Constraints and boundary conditions. Today numerical relativists perform full 3+1 
dimensional simulations with some success. One would think that simulating a static 
problem, namely ^^gravito- statics" , would be well-understood by now, but this turned out 
not to have been the case. The main conceptual hurdle which was necessary to cross 
was the treatment of constraints and boundary conditions. This problem was solved by 
Wiseman in the 2d case [57|, as we now describe. 



Relaxation and electro-statics. Since Newtonian gravito-statics is equivalent to electro- 
statics it is useful to recall the method there. In electro-statics one wishes to find the 
electro-static potential $(x) defined over some domain, satisfying 

• The Poisson equation 

A$ = -47rp, (4.3) 
where p is the given charge density distribution 

• Boundary conditions: Dirichlet, Neumann or some mix. 

A successful numerical algorithm to solve this problem is the relaxation method. This 
method is very physical in the sense that it has some similarities with the way in which 
an excited field settles down or "relaxes" as a function of time to a static solution. In the 
relaxation method one chooses a grid, consisting of points Xij, and then one starts with 
an initial field configuration ^f'j = ^^^\xij) which satisfies the boundary conditions. At 
each step $ is modified according to a local rule to create a sequence which converges 
to the solution as m — > oo. More specifically, once chooses a discretization of the Laplacian 
and solve for ^ij. For example, if one uses a square grid with spacing h and chooses the 
following discretization 

(A^)disc = ^ m - ^] 



^® Attributing {d — 2)(d — l)/2 — 1 degrees of freedom to gravitons in d dimensions. 
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then after discretizing ( [4. 3D using (|4.4| ) and solving for j we find 



(m+l) 



(4.5) 



Actually, one can go further and introduce a "relaxation speed parameter", cu, defined 



by 



^(m+l) ^ ^(m) + ^ 



(4.6) 



such that Lo = 1 corresponds to the rule ([4.5|), while u; > 1 is called "over-relaxation" and 
w < 1 is called "under-relaxation" . Clearly, the solution is a fixed point of the process, 
irrespective of the value of lo. Its importance lies in changing the convergence properties. 
For some interval of oj containing iv = 1 the process is guaranteed to converge since at each 
step the energy is reduced and the solution is a unique and global minimum of the energy. 
In this range uj may be adjusted for convergence speed. 

Gravito- statics. Similarly to electro-statics, General Relativity allows for relaxation. In 
our case there are 5 equations of motion, and after fixing the gauge the equations are split 
to 3 equations of motion and 2 constraints ( "gauge fixing constraints" ) . A convenient and 
quite natural gauge choice is the "conformal" gauge 



The action and constraints in this gauge are as follows. The action is given by 



levr Gd 



where 



dr dz e 



1- 



K^p di<pf' + {D-3){D-A)e 



2{B-C) 



(4.7) 



(4. 



K. 



a/3 : = 



A + {D-2,)C 
A 
B 

C ^ 

1 D-3 

1 D-3 
D-3 D-3 (D-3)(D-4) 



(4.9) 



To express the constraints compactly it is convenient to define 

{Cnsr)ij := didj"^ + dj"^ - K^p dj^f^ 
in terms of which the constraints are 



(4.10) 



{a^Y' {Cnsr)ij 
{a^f {Cnsr)ij , 



(4.11) 
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where a^'^ are the usual PauH u-matrices 



1 



1 



.3 . 



1 



(4.12) 



a' 



1 



-1 



In this gauge Wiseman |57] was able to formulate gravito-statics as a relaxation prob- 



lem. The 3 equations of motion are elliptic of the form 



AX = Src 



(4.13) 



where X is any of the fields A, B,C, A = dzz + drr is the flat space Laplacian, and Src are 
some non-linear functions of the fields. As such we would like to subject them to relaxation 
with some b.c. at infinity and at the horizon. Normally 3 elliptic equations require three 
boundary functions as data. However, the analogy with electro-statics, where one needs 
to specify only the electro-static potential on the boundary, leads us to expect the b.c. to 
consist of a single function. So the question is how to determine the b.c. Another problem 
is how to guarantee that the remaining two constraints, which are hyperbolic in this case, 
are satisfied as well. 

The problems of b.c. and constraints are solved simultaneously by noting the con- 
straints' Bianchi identities. In this specific ansatz the constraints are Grz and Grr — Gzz, 
where G is the Einstein tensor. The Bianchi identities are 



where g = det g^jy and one uses the Einstein tensor with mixed upper and lower indices. 
These are Cauchy-Riemann equations and hence 



is analytic in the r + i z variable. For an analytic function to vanish in a domain, we know 
that it is exactly enough to impose that its real part vanishes on the boundary and that 
its imaginary part vanishes at a point. Alternatively, one may choose and arbitrary a and 
impose the boundary constraint Re(e*" G) = together with Im(e*" G) = at a single 
boundary point. It turns out that these boundary conditions chosen by the constraints' 
Bianchi identities, do not only guarantee the constraints but are also exactly sufficient for 
the elliptic problem. Note that altogether the b.c. consist of a single function, as expected 
from the analogy with electro-statics (plus a function at one point), and that there is some 
freedom in specifying the b.c. (such as choosing a) which is analogous to the choice of 
Dirichlet/ Neumann. 

Unlike electro-statics the action is neither positive-definite nor is it quadratic in the 
fields and therefore convergence is not guaranteed. In practice, however, this method 
performs quite well (see also the "convergence" part in subsection ^^ ). 




(4.14) 




(4.15) 
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Since the Cauchy-Riemann equations are special to 2d it is not obvious how to gener- 
ahze this structure to higher dimensional gravito-statics, which is required in order to solve 
for black holes in backgrounds with higher dimensional compact manifolds {p > 1). One 
may speculate though, that the general rule is the same as the last emphasized sentence. 

Another issue of boundary conditions, which is dealt with in the conformal gauge 
( |4.7| ), is fixing the boundary. In GR, due to the dynamic nature of the metric one cannot 
anticipate the final location of the horizon. However, it is normally best for numerics to fix 
the location of the horizon in coordinate space. The conformal gauge allows for analytic 
coordinate change as residual gauge. It was shown that this freedom exactly suffices in 
order to fix the horizon to be a circle in the (r, z) plane, by writing down elliptic equations 
for the coordinate transformation [^8|. The circle's radius, ph, is a free parameter that 
determines the size of the black hole ( "the holomorphic invariant of the domain" ) . 

Harmark-Obers coordinates 

By choosing the coordinates with some physical logic Harmark and Obers found an 
ansatz with 2 rather than 3 fields p^. The ansatz was constructed for black hole solutions. 



it explicitly fits the uniform string, in |58] it was used for the perturbative analysis around 
the GL point, but it was not used in a relaxation algorithm so far. 

[ |6^ ] start by defining orthogonal coordinates {R, v) over the (r, z) plane such that 
they interpolate between spherical coordinates near the origin and cylindrical coordinates 
asymptotically. To that purpose R was defined through $Ar(r, z), the Newtonian potential 
of a point source at the origin 

oo 

'^n:=p{^-' Y1 (4-16) 

n=-oo (r2 + (z + nL)2) 2 

(see figure , by'^^ 

Kq 



« . (4.17) 

K 

V is defined to be an orthogonal coordinate of period 27r and it can be parameterized such 
that near the origin is approaches x (see figure ^) while asympototically it approaches z. 
This is achieved by the system 

d^v = -kr'^-^dr^N , (4.18) 

where k is some constant and the system is integrable since <I>7v is harmonic. We see here 
that this construction is special to 2d, and it is not clear whether it can generalized to 
higher dimensions. 

The {R, v) plane is drawn in figure |l^. It is a semi-infinite cylinder with one marked 
point (denoted by x), which is (r, z) = (0, L/2), where the equipotential surfaces turn from 



^"jGTj find an explicit expression for $]v by separating variables and expressing the radial functions in 
terms of modified Bessel functions of the second kind. 
"^jeTl normahze R such that i? ~ 27r r/L for r > L. 



-38- 



Figure 18: The equipotential lines of tlie Newtonian potential around a point-like source in 5d 
drawn in the (r, z) plane. Note the change in their topology from spheres to cylinders. 
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Figure 19: The Harmark-Obers coordinates R,v i? is a function of the Newtonian potential 
and V is periodic with period 2tt. The x marks the point {R, v) = {Ri, tt) which is a singular point 
of the coordinate transformation. For R < Ri equal- i? surfaces are spheres while for R < Ri they 
are cylinders. The size of the black hole is changed by adjusting i?o- 



spheres to cylinders and the differential of the transformation from (r, z) to {R, v) vanishes. 
In these coordinates Ro, the location of the horizon, is a free parameter that determines 
the size of the black hole, while the location of the marked point remains fixed at some 
Ri. As the black hole grows Rq reaches Ri and from then on the same equations should 
generate string solutions rather than BHs (for some appropriate boundary conditions). 

In the {R,v) coordinates it turns out that two fields A, K suffice (rather than three). 
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according to the following ansatz 
ds^j^^ = f dt^ 



f-^AdR^ + 



A 



dv^ + KR^ dnl_2 



where 



/ = 1 



r: 



d~3 



Rd~3 



(4.19) 



(4.20) 



The justification for this ansatz was demonstrated (in 6d) |69[] by showing that for any 
solution the equations for a change of coordinates to the ansatz ( 4.19| ) are self-consistent. 
The proof was then generalized in to any dimension. Another way to see that reduction 
is to write a fully general ansatz by adding a third field C, say as the coefficient of dv"^, 
and then one finds that the Gu constraint gives an algebraic relation between the three 
fields |7C]. Yet, it is not clear how the reduction of fields could have been anticipated and 
whether there are other ansatzs with this same property. 

A further reduction is possible, expressing A in terms of K and its second derivatives 
(eq (6.6) of ||6^). However, this step makes the equations of order higher than second, 
and hence unamenable to relaxation (actually even the equations for both A, K were not 
brought to relaxation form so far). 

4.2 Numerical issues 

We turn to a brief overview of the numerical implementation issues of our gravito-statics 
problem. The first simulation of the system was |57] which found the correct prescription 
for b.c. and constraints and applied it to obtain solutions for non-uniform strings including 
highly non-uniform ones very close to the merger point. That work benefitted from insights 



gained from works on the "black hole on a brane" problem such as [71|. A simulated 
black hole on a brane appeared in [^] while the first simulation for caged BHs appeared 
simultaneously in |3^] (5d) and ||3^ (6d). However, due to convergence problems large 
black holes were not possible to obtain. These problems were largely overcome in (5d 
&: 6d) whose figures summarize the state of the art. 



Implementation. The most important decisions for implementing the numerics are the 
choice of coordinates, fields and grid. So far the coordinates and fields were chosen in 



accordance with the conformal gauge (4.7), with variations intended to extract the singular 
and asymptotic behavior for smoother numerics. For example, the ansatz used in |57| for 
6d strings was 



ds' 



G^m + r' 



e^ ^ dt^ + ^ (dr^ + dz^)+e^^ {G^m + 



'dn^ 



(4.21) 



In 1 32 1 the implementation proved to be very sensitive to a redefinition of fields, for some 
unclear reason, the main obstacle being a simple linear redefinition of fields of the form B' = 
B + C, C' = B — C . The more efficient Harmark-Obers coordinates were not implemented 
so far and it would be interesting to do so. 

When choosing a grid one should consider the following factors: the number of nearest 
neighbors for each point, the density of points and multi-grid issues. Hexagonal lattices are 
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particulary suited for the solution of the Laplace equation in 2d, because one can first use 
some relaxation method just to distribute the grid points according to a prescribed grid 
density (to allow for variable grid spacing), and then weights for the discretized Laplacian 
are uniquely determined from the points' location. However, for convenience all simulations 
so far involved square grids and the density of points was sometimes adjusted by using the 
mapping from grid space to coordinate space. For instance in two grids were used: 
in the near zone the grid was evenly spaced in p and cos(x) (see figure ^ for coordinate 
definition) while in the asymptotic zone it was evenly spaced in log(r), z. The need for two 
grids came from the need for different densities of points in the two regions. The price to 



pay is in the complexities of grid matching. Another grid tool employed in [32| was the 
"multi-grid" where the problem was relaxed successively on finer and finer grids, in order 
to accelerate convergence. 

Another, less critical choice, is the choice of discretization of the differential operators 
(such as (0)). 



Convergence. Having implemented the numerical scheme, one typically presses "enter" 
and prays for convergence. To achieve it one is free to use the "convergence speed" pa- 
rameter (jj (4.6). Often this does not suffice and we enter the domain of black magic. 



One difficult problem that arose from an innocent redefinition of the fields was already 
described. 

Another issue involves non-convergence as a result of problems near the exposed r = 
axis (so it does not apply to strings). This problem is known to be solved effectively by the 



following "Wiseman trick" |71]. The field B appears in source terms which are sensitive to 
errors near the exposed axis and generate instabilities. To cure that one replaces B hy & 
numerically distinct variant, B2, which is meant to be more accurate than B near the axis 
and is gotten by integrating a constraint from the axis outward. 

More generally, one expects to have a relation between physical stability of a solution 
and convergence. However, the method converges nicely for non-uniform strings which 
are believed to be unstable. The convergence is the result of boundary conditions at the 
horizon which force the solution to be non-uniform and which perhaps can be described as 
exerting some "pressure" which stabilizes the string also in the physical sense. It would be 
nice to understand this issue better. 

Tests. The results of the simulations must be tested. The following tests were employed 

• Convergence. 

• Constraints. 

• Integrated first law. 

This law (^3) relates m, r which are measured at infinity with /?, S measured at the 
horizon. As such it provides a strong overall test of the numerics, and it is satisfying 
to find that it holds. However, note that the first law does not test the constraints 
appendix C]. 
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Comparison with analytic results. 



Some analytic results for small black holes (see subsection 4.4) allow testing of the 
numerical data. The satisfying agreement is demonstrated in figures 2C, 21, 2S 
taken from |32|. In all figures the size of the BH is measured by x := 2 ph/L where ph 
is the horizon radius in conformal coordinates and this "numerical" definition of x, 
valid in this subsection, should not be confused with a different "analytic" definition 
of the small parameters x valid in subsection 4.4. Stars, diamonds or triangles repre- 
sent actual simulation results while the curve represent smooth interpolations which 
coincide with the theoretical predictions of |25] within the simulation precision. 



4.3 Time evolution 



As remarked in subsection 2.5, "issues", the time evolution of the unstable system raises 
deep questions regarding singularities. Thus there is ample motivation for a dynamical 
numerical simulation. On the other hand, we expect any such simulation to crash before 
the singularity is reached. So in order to achieve progress it is probably necessary to have 
a theoretical local model for the time evolution, which could then be tested and supported 
by a numerical simulation. 



Thus far a single time evolution simulation was performed: |72| simulated the decay 



of a sub-GL uniform string in 5d. The main features of the evolving spacetime can be read 



from figure 24 . The initial configuration is an unstable string of radius rg on top of which 
the unstable perturbation grows in time. The region where the horizon moves radially 
inward (the "waist") collapses fast, stretching at the same time in the z direction (growing 
gzz) until the minimal (areal) radius reaches about 0.15 ro when the grid stretching is so 
large that the simulation cannot proceed. In the region where the horizon grows (the "hip") 
the metric approaches the metric of a 5d BH with a comparable mass, so that the maximal 
(areal) radius at the horizon is about 2ro. 

Opinions on the interpretation of these results differed. The author of this review 
interprets them as strong evidence against the formation of a stable non-uniform string 
end-state. The authors of [^2| are much more careful and would only say that the results 



- 42 - 



0.08 r 
0.06 
i 0.04 
0.02 
0- 



(1 = 1.48 x^-1.31 







0.1 



0.2 



[> 



0.1 



,10" 



0.2 



6 
4 
2 

^ ^0 



Q 3 Ag = 1 9.54 x^+5.6 x'* 



0.2 



0.1 




0.05 0.1 0.15 0.2 0.25 

X 



1/T=6.25x + 2x'' 



V 




0.1 0.2 

X 



i2 
1.5 
1 

0.5 

-0 



Figure 21: Thermodynamic quantities, ^, t, ^3, 1/T which arc the mass, tension, area and inverse 
temperature, respectively, of a 5d black hole as a function of its size. Simulations and perturbative 
analytic expressions agree, except for the tension results which are less reliable. Reproduced from 
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Figure 22: The area in units of the surface gravity, as a function of 5d black hole size. This 
quantity measures a correction to the temperature due to the non-zero potential at the origin from 
the image BHs. Simulation and perturbative analytic expression agree. Reproduced from ||3^ . 



"are not inconsistent with Gregory and Laflamme's conjecture that the solution bifurcates 
into a sequence of black holes ... [On the other hand it is] not necessarily inconsistent 
with [43|... At the same time, a continuation of the observed trend would argue against 
achieving a stationary solution with a mild dependence on the string dimension". 

Later [^3| presented a re-analysis of the metric obtained in [^] which somewhat clarifies 
the picture. observed that the affine parameter on the horizon grows extremely fast, 
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Figure 23: The inter-polar distance (going around the compact dimension), defined in figure g8[ as 
a function of black hole size in 5d. The line denotes a best-fit curve. Its quartic nature is explained 
in pUA, though a theoretical prediction of the prefactor is not presently available. Reproduced from 




Figure 24: A time sequence of embedding diagrams for the event horizon of a decaying black string 
in 5d, from a numerical time evolution. Initially, the string is nearly uniform (top), and by the 
time the simulation comes to a stop due to grid stretching (bottom) , the horizon becomes that of a 
black hole with a thin pipe connecting the poles (going around the extra dimension). Reproduced 
with permission from ||7^. 



faster than the exponential of asymptotic time, and suggested that the horizon might pinch 
off in infinite affine parameter. This possibility was mentioned in but argued against. 

In light of the critical dimension D* = "13.5" it would be very interesting to run 
the simulation again for a range of dimensions, and to test/ confirm that for D > 14 the 
unstable string settles down quickly to a slightly non-uniform string. 

Method. Let us count the number of degrees of freedom, which is of course larger 
than in the static case. The domain is the 3d space parameterized by (r, z, t) and the fields 
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are the metric and the scalar which is roughly qqq or C in (^Jj). The total number of 
fields is 7. After choosing some gauge in the 3d domain the number of fields reduces to 4 
(compared to 3 in the static case). Finally, the "physical" number of degrees of freedom is 
a total of + 1 = 1 (just like the static case): none for the graviton and +1 for the scalar. 
The most general time-dependent metric is 

ds^ = (-a^ + jABf3'^f3^)dt^ + 2^AB[i^dx'^dt + -lAsdx'^dx^ + -f^d^'^ , (4.22) 

where = (r, z) , and d^l"^ is the 2-spherical line element with coordinates chosen orthog- 
onal to the t = constant congruences (hence there is no shift corresponding to angular 
directions). All metric components in ( [4.22| ) are defined in the 3d domain. [|72| choose the 
gauge 

a = ast 

999 = , (4.23) 

where the string metric function is read from the ingoing Eddington-Finkelstein form 

dslt = -(1 - 2M/r) dp + 4 Af/r dr di + (1 + 2M/r) dr"^ + dz^ + d^^ . (4.24) 

They also comment that "In a preliminary version of our code, we also required that 
j3'^ = 13^Q^. This, however, caused a coordinate pathology to develop at late times during 
the evolution of unstable strings — specifically, some regions of the horizons approached a 
zero coordinate-radius, while maintaining finite proper radius." Thus the condition on (3^ 
was replaced by one on ggg. 

Another hurdle turned out to be the boundary conditions at infinity. To eliminate 
cut-off problems spatial infinity was brought to a finite point by replacing r by r/(r + 1). 

4.4 Analytic perturbation method 

While we do not know to write the black hole metric in closed form, metrics for small black 
holes can be well- approximated everywhere: for p L the D dimensional black hole is 
a good approximation, while for p ^ po (po is the Schwarzschild radius) the Newtonian 
approximation is good. Moreover the two approximations have an arbitrarily large overlap 
in the small black hole limit. Therefore, one expects that the black hole metric can be 
systematically expanded in a perturbation series with a small parameter being x := pq/L.^'^ 
Such a perturbative method was studied in several papers: the general procedure and first 
order results were given in ||2^ (pre-announced in [^l|), the full second order in 5d was 



obtained in [29| and for a general dimension the Post-Newtonian order was performed in 



|3C]. A different, one zone approximation was given in [27|. 
The objectives here are to 

• Provide tests for numerics. 



^■^This "analytic" definition of x, valid in this subsection, should not be confused with a different "nu- 
merical" definition of the small parameters x valid in subsection [4. 4 
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Figure 25: The metric for small black holes can be obtained through a "dialogue of multipoles" 
matched asymptotic expansion. The two zones are the asymptotic zone p S> po and the near zone 
p <C i. The smaller the black hole the larger is the overlap region between the two zones. 



• Extrapolate to the phase transition region. 



Method. It is not possible to use here the "usual" perturbation method, the one where 
a "zeroth order" solution is deformed order by order to follow the deformation of a small 
parameter of the problem, since here the domain of coordinates changes with x, and since 
we do not have a zeroth order solution. However, one can use the well-known technique 
of "matched asymptotic expansion". In [28| this technique was applied by defining two 
zones: an "asymptotic zone" and a "near horizon zone" (see figure [25|) , which have a large 
overlap in the limit x ^ 0. The metric is solved perturbatively in each zone, with boundary 
conditions coming from matching with the other zone. The need for matching produces 
an intricate (and dimension dependent) pattern of crossings between the various orders in 
the two zones - "the perturbation ladder" - see figure 26. Effectively the gravitational 
field produced by the images changes the shape of the BH, or its mass multipoles, and 
that in turn back-reacts and changes the field multipoles. This procedure is much simpler 
than the usual dynamic matched asymptotic expansion on account of being static, and was 
therefore given a special name "a dialogue of multipoles" [28|. 

The zeroth order. In the asymptotic zone the zeroth order metric is flat space with a 
compact dimension and the origin removed, while in the near zone it is the Schwarzschild 
black hole, with the periodicity L being far away and invisible. 

The Newtonian approximation. The first order analysis in the asymptotic zone is easy 
- it is just the Newtonian approximation. Here one considers the metric in the weak field 
limit 

g^iv = Tlfiu + h^t, , (4.25) 

where r/^,y is the flat space (Minkowski) metric and hfj_i, ^ 1 is a small correction. One 
defines ^ 

hfiiy ■■= h^y - -h^ ri^y , (4.26) 
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Figure 26: The "perturbation ladder" for tlie "dialogue of multipoles" matctied asymptotic ex- 
pansion. The upper row depicts the asymptotic zone and the lower row the near zone. Each box 
denotes a certain order in the perturbation series - in the asymptotic zone one counts orders of po 
and in the near zone orders of L^^. The arrows denote the flow of information between the two 
zones, in supplying boundary conditions through matching. The figure shows the general pattern 
for an arbitrary spacetime dimension D (d in the figure is denoted by D in this review). Reproduced 
from [|8). 

in terms of which the harmonic gauge^^ is chosen 



df'h^, = . 



(4.27) 



The linearized field equations read 

1 



Dh 



(4.28) 



where □ is the flat space D'alambertian. Solving hfj^u for weak and slowly moving sources 
where the only non-zero component of T^j,y is taken to be Too, and inverting ( |4.26 ) through 



h — h 



(4.29) 



we find 



htt 



N, 



D-3 



(4.30) 



^Also known as the De-Donder or Lorentz gauge. 
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where the Latin indices stand for the spatial components and is a the Newtonian 
potential (|4.16| ).^^ In the D = 5 case there is a useful way to perform the summation in 



the definition of <I>Ar ( [4.16| ), yielding 

- 17 cosh (if^)- cos (^) • ^'-''^ 

Linear perturbations around Schwarzschild . The first correction to the near zone 
requires more work than the Newtonian approximation in the asymptotic zone - it is 



a generalization to higher D of the well-known paper by Regge and Wheeler |74]. The 
metric is given in terms of a single function E = E{p, x) which determines the metric. 
After separation of variables, Ei{p) satisfies a "master equation" 

D'^E f2 1 I \ dE , , X + (D-\)w , , 



\X X-1 X-wJ dX X {X -1){X -w 

where I is the angular momentum and 



\D-3 



D-2 



wld ■-- 



{l-l){l + D-2) 



P,o := ^ . (4.33) 

There are four regular singularities in the complex plane including oo,'^^ so it is a "Heun 
equation", by definition. The singularities are at 0, pojOo and at w, but the latter can 
be eliminated by a field re-definition (re-defining E) and in that sense is "non-physical". 
Since there are 3 remaining singularities the solutions can be expressed in terms of hyper- 



geometric functions. These results were first obtained in [76| and then in |28| each using 
somewhat different methods. 

Regulation of divergences at Post-Newtonian order. At the next, Post-Newtonian order 
divergences were identified and regulated in |^0[ . The divergences originate from integration 
over the non-compact overlap region (which is the neighborhood of the origin for the 
asymptotic zone). The regularization can be defined through a "cut-off and match" method 
where one places a cut-off on the integral, matches with the other zone, and then sends 
the cut-off away to remain with a finite solution. This is equivalent to a regularization 
known as Hadamard's partie finie, and it is closely related to the concept of "subtracting 
the self-energy". This regularization allowed to obtain the thermodynamic quantities at 
Post-Newtonian order. 



^''Possibly up to a constant pre-factor, depending on conventions. 

^^Recall that xq is a "regular singularity" of the differential equation a{x) y" + P{x) y' + 7(a;) j/ = when 
after normalizing the equation by overall multiplication such that a(a;o) 7^ 0, oo either l3{x) or 7(2;) has a 
pole at a::o, but the order of the /3-pole is at most one and the order of the pole in 7 is at most two. 

^''For a short discussion of the Heun equation see appendix A in psl]. More information can be found in 
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In [^] an alternative method was given employing a single patch, and using the effi- 
ciency of Harmark-Obers coordinates |6^]. There a "first order" approximation is given, ^"^ 
and I expect that the method could be developed to a full perturbation series by succes- 
sively improving a suitably chosen initial guess. It has the advantage of using only a single 
zone, and doing away with matching. However, the method depends on the choice of ini- 
tial guess, and unlike the previous method, the differential operator to be inverted would 
probably change at each order. 



Results. In 5d the second order metric was obtained in [29| from which the following 
thermodynamics were deduced 
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where the small parameter is defined through 6 = 8^5 M / (Sir L'^) = Pq/L?'. 

For arbitrary D the best available determinations of the thermodynamics quantities 
|2|, H] are 
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where 



D-3 



(4.35) 



(4.36) 



and C{s) := 1/n* is Riemann's zeta function, which appears here due to the sum over 

black hole images. We note that the black hole tension ( 4.35| ) vanishes to lowest order, 
namely 0(M) and the leading order result can be explained by the Newtonian potential 



between the black hole and its images |30|. 

When one considers results beyond the leading order, such as some of the results 
above, their "scheme dependence" should be borne in mind. "Scheme dependence" is used 
here to mean the freedom to re-parameterize po and through it the branch of small black 
hole solutions: while to leading order all definitions of po coincide (up to a multiplicative 
constant), there is no unique or natural definition for the subleading corrections, and 



''^The first order in tlie metliod of 
multipoles |30[. 



actually includes also the second order results of the dialogue of 
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Figure 27: The first non-spherical deformation of the horizon is a quadruple moment, making 
the black hole prolate, namely elongated along the z axis. In order to measure the eccentricity we 
define the area of the equator sphere and Ax_ the area of a "polar" sphere. 



therefore results should be accompanied by a specification of the scheme in which they 



were obtained. In particular, the definition used in is specified in |28|, subsection 3.3.2. 

Some nice geometric quantities can be measured as well. The leading deviation from 
spherical shape is such that the BH becomes elongated along the z axis, namely prolate - 



see figure 27. Its eccentricity can be measured by 



^± , (ZJ-3)^r^(2 + ^)Cp-l) ^po^D-i 

^ " 8p-2)r(^) (t) • (^-^'^ 

This result actually comes from beyond first order. The inter-polar distance (going around 
the compact dimension, see figure |2l) is given by 



Vies = ^-2^=^/50 + 



r- r(gE|) 



= 4— ^ . (4.38) 

where the ellipsis is order o [{pq/ L)^^"^^ ■ L. Since Is = in 5d the black hole makes 
room for itself, exactly compensating its size, and for that reason it was called "the black 
hole Archimedes effect" |3^]. More precisely, in 5d L^^^^^/L = o (^{po/L)^), and numerical 
results indicate that indeed the next order is non-zero: L-^^^^^/L = O [{po/L)^) . For higher 
D, < Id < 1, the effect is milder, and in addition, the order of the ellipsis depends on a 
choice of scheme defined by eq. (B.15) of 



Impact on objectives. The ability to produce results for small black holes from two 
quite different methods, perturbation theory and numerical analysis, allowed for a healthy 
feedback between the two, testing and perfecting each one of them - see subsection W3 and 



especially figures EG, 21, 22, 031 While the perturbative method does not apply directly to 
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Figure 28: The inter-polar distance is defined to be the proper distance between the poles, going 
around the compact dimension. The figure depicts the tendency of mass to "make room for itself" , 
an phenomenon termed "the black hole Archimedes effect" in [p2|. 



the phase transition region, it provided important tests for the numerics which are capable 
to study that region. 

Implications for the phase diagram. Even though the phase transition region is outside 
the domain of validity of the perturbative method, the results above can be extrapolated to 
provide evidence [^0| for the critical dimension separating first and second order behavior, 
as we now proceed to explain (of course chronologically D* was discovered first by [^] ). 
By extrapolating the result for the tension of the small black hole ( [4.35| ) one obtains an 
extrapolated straight line in the phase diagram whose intersection with the 6 = axis, 
called should be compared with ^gl the critical Gregory-Laflamme mass. While for 
small D numerically one finds that > fJ-GL, for large D examination of the asymptotic 
expansions show that fix ^ l^'GL, which means that for large D the black hole phase is 
not available yet at the onset of the GL instability and therefore the string must decay 
to a different phase, which is evidence for the smooth second order decay into the stable 
non-uniform string. 



5. Related work 



Here we mention some related work. 



Relation with gauge theory & charged black holes. The gravitational phase tran- 
sition is closely related to a gauge theory phase transition, the so-called "Gross- Witten" 
transition where the eigenvalues of some unitary matrix (originally a Wilson line over 
a plaquette, in a later application ^8|, a Wilson line around a compact dimension) change 



from a clumped distribution to a uniform one as the temperature is raised (see figure 29). 
This correspondence was first noted by Susskind (1997) ||7^ on a somewhat qualitative level 
and was fully analyzed quantitatively by Aharony-Marsano-Minwalla- Wiseman in . The 
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Figure 29: In the Gross- Witten phase transition [[77| the dynamics of the eigenvalues of a unitary 
matrix in gauge theory is mapped to a system of particles in a horizontal loop in the field of gravity. 
In large N (the number of particles) the particle density becomes a continuous function p{9). The 
phase transition is between localized to non- uniform and then to uniform distributions. Susskind 
(1997) identified this phenomenon as the gauge theory dual of the black-hole black-string transition. 



correspondence with gauge theory enriches the system by adding another parameter to the 
problem - the coupling. 

In |^0| the authors showed that the phase diagram of black holes/strings directly maps 
onto the phase diagram of near-extremal black branes, and hence via holography onto the 
phase diagram of the dual non-gravitational theories (Super- Yang-Mills and "Little String 
Theory" ) . More specifically, they present the phases of the non-gravitational theories and 
analytical results for the thermodynamics of the localized phase and the new non-uniform 
phase (connected to the "image" of the GL point). In fact part of the early motivation 
of these authors for studying black holes in Kaluza-Klein backgrounds (see ||6^) has been 
precisely that, namely to relate it to gauge theories. See also the brief review pT| . 

In 1 82 1 the GL analysis was generalized to "point-like" charged strings^*^, and the 



instability was found to persist for all charges, and moreover kci diverges as extremality 
is approached. 

For ideas about the relation of this phase transition to the Hagedorn phase transition 



in string theory see |83| and references therein. 



Bubbles. Taking the metric of the black string (2^) and performing a double analytic 
continuation in z,t (namely t' = iz, z' = it) one gets an unstable metric called the "bubble". 
In the Euclidean setting the analog would be to exchange /3 <-> L thereby inverting 
which is not that interesting. From the Lorentzian perspective this is another phase that 
could be put on the same phase diagram. It has the special property t L = {D — 3)M. 
Thus the bubble exactly saturates the Strong Energy Condition upper bound on tension 
(see subsection |3.1D , which motivated [51| to study this phase as well as combinations of 



several bubbles and BHs. So far bubble phases were not seen to connect with the phases 
discussed in this review, so it is conceivable that the two issues are decoupled. See also the 



brief review [Sl|. 



Braneworld black holes. The topic of black holes on braneworlds is closely related to 
the subject of this review. There the dimensionless parameter is the ratio of the black hole 

**In 5d a string may carry either point-like or string-like charge under a Maxwell field. 
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size and the size of AdS. A comprehensive bibliography is beyond the scope of this paper, 



but we shall only mention braneworld work |71, 84, 
compact dimensions. 



85 1 which was followed by work on 



6. Summary and Open questions 

Why review. We start with a description of some of the major recent developments in 
the field from a personal perspective, leading to the decision that the field is ripe for this 
review to be written. 

The Gregory-Laflamme instability and the whole black-hole black-string transition 
have been attracting constant attention in the String Theory community as well as in 
General Relativity since its discovery (1993).^^ In may 2001 the interest rose further due 



to the provoking ideas in ||4^. A year later I suggested in [44| a qualitative form for the 
phase diagram which proposed the merger of the black hole and black string phases, and 
rested on Morse theory and a topology change analysis. Even today that paper includes 
most of what I know about the system, but back then many felt that the evidence was 
not convincing enough, and the paper was even refused for publication by two leading 
journals. The non-trivial nature of these predictions is best illustrated by the various other 
possibilities that were considered in the literature, see for example the six scenarios in p6| , 
section 6. 

Now,^'^ after more than two years and more than 30 papers, we have a much more 
detailed knowledge of the phase diagram, culminating in the numeric results of |^4| which 
combine new numerical data for black holes with older data for non-uniform strings into 
a full phase diagram (figure 16). This full, actual-data phase diagram exhibits the merger 
of the phases, vindicating concluding most of the uncertainty around |43, 44] as we 
discuss further below, and providing a natural point in time to write a review. 

Below we summarize the results and open questions. The results include two big 
surprises: critical dimensions and topology change, but we remind the reader that while 
there was much progress in understanding the static phase diagram, there was practically 
no progress towards the deeper issues regarding time evolution and singularities. 

6.1 Results 

• Confirmation of emerging phase diagram. 

The current understanding of the qualitative features of the phase diagram^^ is sum- 



marized in figure 15, where the D < 13 case originates in H]. We note that all 



diagrams include a merger point, 
''^with a steady flux of about 20 citations/year. 

^"Refers to November 2004 when the first version of this review appeared in the archives. 
Again, we stress that this is not actual data, but rather what seems to be the simplest possibility given 
the available knowledge, some of which was obtained through research designed to obtain these qualitative 
features. 
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For D < 13 the end-point of decay is seen to be a black hole, and not a stable black 
string which was predicted in I^Sj. The evidence for that includes a numerical de- 
termination of practically the whole branch of the non-uniform string which emerges 
from the GL point and was found to have higher mass than the critical string and 
thus cannot serve as an end-point for decay |57]; the continued collapse of a time 
evolution simulation coming to a stop due to grid stretching, rather than the 
predicted stabilization; theoretical difficulties in proposing a phase diagram which 
includes the predicted phase and satisfy the Morse theory constraint 1 44 ] ; and finally 
the fact that such a phase was not found neither analytically (see the attempt of [^] ) 
nor numerically. 

For D > 14 the end-point of (smooth) decay does turn out to be a stable non- 
uniform string (as a result of |^^), thus partially vindicating However, I claim 
that in a deeper sense, this case does not validate the arguments of [^]. First, the 
arguments were independent of dimension, and second, even for D > 14, horizon 
pinching, which finds to be forbidden, must happen after the non-uniform string 
evaporates enough. 



Thus while the claim of |4S] stimulated much of the research reported here, and 
while there is no doubt regarding the calculations presented there, we find strong 
(actually, overwhelming in my opinion) evidence against it. One should therefore 
reconcile the derivation in with the results presented above. However, it is clear 
that if something went wrong in the arguments of [^] it is rather deep and there are 
lessons to be learned from it. Following the analysis of the numerical time evolution 
in [^] it seems that while indeed pinching does not occur in finite horizon affine 
parameter, consistent with the findings of |43|, it could well happen with infinite 
horizon affine parameter and finite asymptotic time (see also ^^). This possibility 
was indeed discussed in |43], and finally rejected, but with less certainty. In 



this point of view seems to have been accepted: "It was conjectured by Kol [14] 
that the nonuniform black strings should meet the squashed black holes at a point 
corresponding to a static solution with a singular horizon, and this appears to be the 

case [15]. ...the black string horizon might pinch off in infinite affine parameter 

[the] suggestion ... that the black string will break up into spherical black holes might 
still be correct." Altogether, it looks like a reasonable consensus arises where the 
central calculation of |4^ holds and only the arguments against pinching in infinite 
time had some loop-hole, which invalidates, however, their bottom-line conclusion 
regarding the end-point of the decay into stable non-uniform strings. 

Another possibility for a flaw in the assumptions of [43| should be mentioned, namely 
the reliance on the increasing area of the event horizon. Even though the authors 
were careful to use a version of the theorem valid even if there is a singularity on 
the horizon, it is still quite possible that in the present context the whole notion 
of the event horizon is ill-defined and/or that the singularity completely leaves the 
horizon. The problem with the definition of the event horizon is that if a singular 
shock wave emerges from a naked singularity and reaches null infinity, then we do 
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not have the asymptotically Minkowski null infinity which is normally used in the 
definition. Clearly, this issue, like others involving the time evolution, is not well 
understood yet. 

Critical dimensions. 

The stability of cones, which have a role in the local model for merger, exhibits a 



critical dimension -Dmerger = 10 IH' see subsection |3.4 At a different point in 
the phase diagram, the GL instability changes from explosive (first order) to smooth 
(second order) for D > Dq^ = 13.5, as was shown by Sorkin For the canonical 
ensemble this happens at D > D^j^ ^ = 12.5 |Q. 

These critical dimensions are fascinating, but it is not known whether they have a 
deeper, wider meaning. The critical dimension of the BKL theory of the approach 
to a space-like singularity, D* = 10, is closely related. On the one hand one cannot 
help wondering whether there could be a connection with the critical dimension for 
super-strings, lOd, but on the other hand the derivations and issues involved differ 
widely, for example, in the phase transition no supersymmetry is involved while it is 
central for the super-string. 



The "merger" topology change |44]. 



The spaces of metrics for two different topologies, the Euclidean black string and 
the black hole (in D > 5) are glued together (subsection |3.4|) and a line of Ricci flat 



solutions connects them, see subsection 3.4 (in finite distance, see appendix 



Additional results include 

• The use of Morse theory in GR [44], described in subsection 3.3. 



The role of tension and the first law for this system [46, 31], subsection 3.1 



Formulating Gravito-statics by relaxation (in 2d) [57[, subsection 4.1 



Numerical solutions: strings |5^ and black holes [^, subsection and a 

dynamic time-evolution [^], subsection 4.3 . 

Developing an analytic perturbation method for small caged black holes - "a dia- 



logue of multipoles" 
subsection liT 



29|. See ]27] for a different analytic method. Described in 



6.2 Open questions 

• The deepest questions remain: are there a naked singularity and perhaps, a singular 
shock wave in the time evolution? - see the discussion in subsection 2.3. This problem 
could benefit from both a theoretical approach and from numerical analysis. 



Obtain solutions with discrete self-similarity (DSS) which were conjectured in ]63] to 
describe the merger solution locally near the pinch point (in progress). 
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Gubser and Mitra [^6[ conjectured already in 2000 a connection between perturbative 
and thermodynamic instabilities of strings which is closely related to the issues of this 
review. It found confirmation in several cases, some recent ones being |89, 91| 



(in the last case what appeared to be counter-examples became consistent with the 
conjecture after appropriately generalizing it to allow for scalar fields). It would be 
interesting to settle this conjecture. 

Formulate gravito-statics in more than 2d. 

Numerically trace the phase diagram also for additional dimensions in the ranges 
10 < L> < 13 and D > 14. 

Run a time evolution for D > 14 to observe the second order transition. 



Two items present when the first version of this review appeared in the archives, are 
by now mostly resolved 

• Determine the critical dimension for torus compactifications. 

— torus compactifications do not reduce the critical dimension [^] . 

• Stability analysis in the micro-canonical ensemble. In the canonical ensemble the 
Morse analysis, which identifies the thermodynamic potential F with a Morse func- 
tion, allows us to read the stability of phases off the phase diagram. It is plausible that 
a generalization to other ensembles exists. For instance, that would allow to settle 
the conjectured (micro-canonical) instability of the non-uniform string (for D < 14). 

- Poincare's method (see for a review) teaches us that indeed the stability in 
the micro-canonical ensemble can be determined from the form of the phase diagram 
just like in the canonical ensemble, since stability changes can occur only at turning 
points or at vertices where several phases meet. 



The ring. We list some questions regarding the physics of the sister system of rotating 
rings, even though they are not part of the main subject of this review. 

• Obtain the full phase diagram: stability of phases, order of transition, critical points. 

• Solutions in D > 5? 
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A. Formulae for action manipulation 

Here we collect some formulae which are useful for manipulation of actions. 
The Ricci scalar in the presence of a general fibration 

ds"^ = dsjc + Y^ e^^' dsli 

i 

R = Rx + Yl [^~'^" - 2 ^(^^) - di{dF,f] -Y,di dj (dFi ■ OFj) (A.l) 

i i,j 

where the fibration fields depend only on the x coordinates Fi = Fi{x), Rx, Rvi are the 
Ricci scalars of the spaces X, 1^, di are the dimensions dim(li), and the Laplacian (A) 
and grad-squared {d ■ d) are evaluated in the X space. 

The Ricci scalar of a conformally transformed metric (see for example [p2| ) 

ds^ = e^^ds^ => 

R = e-^'"[R-2{d-l)Aw-{d-l){d-2) (dw)^] (A.2) 

where d is the dimension of the space and the Laplacian and grad-squared are evaluated 
in the non-tilded metric. 



B. Topology change is a finite distance away 



Scaling down a smoothed cone gives a family of metrics which approaches the singular 
cone, as discussed is subsection |3.4| (see figure 0). In this appendix we wish to show that 
the singular cone is at a finite distance in moduli space, just like the conifold. 

Let us denote a smooth cone metric with some specific length scale (of the smoothed 
tip) by 

= dp'^ + e^^'^p^ dnlr. + e^^^p^ dnl„ . (B.l) 

The family of rescaled cone metrics is defined by 



ds' 



-2a 



ds 



to O" = OO (fj 



(B.2) 
OO is the 



and we wish to compute the distance in moduli space from a 
non-smooth cone). 

The metric on moduli space can be found by adding an auxiliary coordinate t, making 
a t-dependent, a = cr(t), and evaluating its kinetic term. Since one wants to hold the 
asymptotic form of the cone fixed it is useful to introduce p = p. Then the rescaled 
cone metric is ds^ = dp^ + g^ dJ^^-he^ (^(/"^")-'^) dQ^ , where dQrn = dnl^. Adding 
i-dependence and substituting back to p we get 



ds^ = dt^ + e""^" {dp- dp dtf ^ e^" dQrn + * 

Using ( A.l| ) the Ricci scalar is 

-m(m + 1) {d{a - a)f - 2m A(a - a) 
-n{n + l){d{b- a)f -2nA{b- a) 
-2mn d{a — a) d{b — a) , 



(B.3) 



R 



(B.4) 
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where the differential operators are in the {p, t) plane. All terms which vanish when d = 
must cancel since the resolved cone is Ricci-flat, and one finds 



R= -2 [{m a + n b')p - {D - 1)] d 



(B.5) 



p^{m{m- l)e-^'' + n{n-l)e-^'') + 2p{D -l){ma' + nb') -D{D-l) 



A useful test for this expression is that it must vanish when substituting the singular cone 
metric rather than the resolved cone. 

The resulting bulk action after multiplying by ^ = e'^"+"^~^'^ and eliminating the 
second derivatives using integration by parts is 



-p^{m{m - 1) e"^'* 



(B.6) 



+ n(n — lie 



2 p{ma' + nh') + D{D - I) 



This action should be regularized by comparison with the singular cone. This can be 
done by introducing a large-distance cut-off p = A and therefore p\ = Ae'^ and subtracting 
for the cone. We are interested in large a and hence p\ ^ oo. In this limit the integral 
seems to be dominated by the large p behavior of the integrand: the large p behavior of 
a, b compared to the singular cone is given by the linearized perturbations ( 3.31| ) 



6a,6br^p', Re{s) = -{D -2)/2 
The large pA behavior of the kinetic term in ( |B.(]| ) is 

rAe" 



'Da 



dppD-^p-i^-^) 



Hence the metric on moduli space in the large a limit is 



and the distance is finite 



exp( — a j da < oo 



(B.7) 

(B.8) 

(B.9) 
(B.IO) 



for D > 2, namely always, since we were only interested in D > 5. 
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